{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Using Command-Line Interface (CLI) from Notebook"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "This example shows how to use the ANDES CLI from Jupyter Notebook.\n",
    "\n",
    "It is based on the iPython magic `!`. To run a shell command from within iPython or notebook, place a `!` before the command."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "For example, to list the directory , use `!ls`."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "1. simulate_and_plot.ipynb       6. using-cli-from-notebook.ipynb\r\n",
      "2. inspect_data.ipynb            kundur_full_out.lst\r\n",
      "3. eigenvalue.ipynb              kundur_full_out.npy\r\n",
      "4. inspect_models.ipynb          kundur_full_out.txt\r\n",
      "5. profiling.ipynb\r\n"
     ]
    }
   ],
   "source": [
    "!ls"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Likewise, to run `andes`, use `!andes`. Addition arguments can be passed as usual."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "ANDES 0.6.9 (Git commit id c8cba088, Python 3.7.6 on Darwin)\r\n",
      "Session: hcui7, 02/23/2020 02:16:10 PM\r\n",
      "This program comes with ABSOLUTELY NO WARRANTY.\r\n",
      "\r\n",
      "usage: andes [-h] [-v {10,20,30,40,50}] {run,plot,misc,prepare,selftest} ...\r\n",
      "\r\n",
      "positional arguments:\r\n",
      "  {run,plot,misc,prepare,selftest}\r\n",
      "                        [run]: run simulation routine; [plot]: plot simulation\r\n",
      "                        results; [prepare]: run the symbolic-to-numeric\r\n",
      "                        preparation; [misc]: miscellaneous functions.\r\n",
      "\r\n",
      "optional arguments:\r\n",
      "  -h, --help            show this help message and exit\r\n",
      "  -v {10,20,30,40,50}, --verbose {10,20,30,40,50}\r\n",
      "                        Program logging level. Available levels are 10-DEBUG,\r\n",
      "                        20-INFO, 30-WARNING, 40-ERROR or 50-CRITICAL. The\r\n",
      "                        default level is 20-INFO.\r\n"
     ]
    }
   ],
   "source": [
    "!andes"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Run a simulation"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "ANDES 0.6.9 (Git commit id c8cba088, Python 3.7.6 on Darwin)\n",
      "Session: hcui7, 02/23/2020 02:16:11 PM\n",
      "This program comes with ABSOLUTELY NO WARRANTY.\n",
      "\n",
      "Parsing input file </Users/hcui7/notebooks/kundur/kundur_full.xlsx>\n",
      "Input file kundur_full.xlsx parsed in 0.3180 second.\n",
      "-> Power flow calculation with Newton Raphson method:\n",
      "Power flow initialized.\n",
      "0: |F(x)| = 14.9283   \n",
      "1: |F(x)| = 3.60859   \n",
      "2: |F(x)| = 0.170093  \n",
      "3: |F(x)| = 0.00203827\n",
      "4: |F(x)| = 3.76414e-07\n",
      "Converged in 5 iterations in 0.0068 second.\n",
      "Report saved to </Users/hcui7/notebooks/examples/kundur_full_out.txt> in 0.0006 second.\n",
      "-> Time Domain Simulation:\n",
      "Initialization tests passed.\n",
      "Initialization successful in 0.0176 second.\n",
      "  0%|                                          | 0/100 [00:00<?, ?%/s]<Toggle 0>: Applying status toggle on Line idx=Line_8\n",
      "100%|███████████████████████████████| 100/100 [00:00<00:00, 152.67%/s]\n",
      "Simulation completed in 0.6550 second.\n",
      "TDS outputs saved in 0.0154 second.\n",
      "-> Single process finished in 1.2587 seconds.\n"
     ]
    }
   ],
   "source": [
    "!andes run ../kundur/kundur_full.xlsx -r tds"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Check the output `lst` file"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "     0,                  Time [s],                         $Time\\ [s]$\r\n",
      "     1,            GENROU delta 0,                 $\\delta\\ GENROU\\ 0$\r\n",
      "     2,            GENROU omega 0,                 $\\omega\\ GENROU\\ 0$\r\n",
      "     3,              GENROU e1q 0,                   $e'_q\\ GENROU\\ 0$\r\n",
      "     4,              GENROU e1d 0,                   $e'_d\\ GENROU\\ 0$\r\n",
      "     5,              GENROU e2d 0,                  $e''_d\\ GENROU\\ 0$\r\n",
      "     6,              GENROU e2q 0,                  $e''_q\\ GENROU\\ 0$\r\n",
      "     7,            GENROU delta 1,                 $\\delta\\ GENROU\\ 1$\r\n",
      "     8,            GENROU omega 1,                 $\\omega\\ GENROU\\ 1$\r\n",
      "     9,              GENROU e1q 1,                   $e'_q\\ GENROU\\ 1$\r\n",
      "    10,              GENROU e1d 1,                   $e'_d\\ GENROU\\ 1$\r\n",
      "    11,              GENROU e2d 1,                  $e''_d\\ GENROU\\ 1$\r\n",
      "    12,              GENROU e2q 1,                  $e''_q\\ GENROU\\ 1$\r\n",
      "    13,            GENROU delta 2,                 $\\delta\\ GENROU\\ 2$\r\n",
      "    14,            GENROU omega 2,                 $\\omega\\ GENROU\\ 2$\r\n",
      "    15,              GENROU e1q 2,                   $e'_q\\ GENROU\\ 2$\r\n",
      "    16,              GENROU e1d 2,                   $e'_d\\ GENROU\\ 2$\r\n",
      "    17,              GENROU e2d 2,                  $e''_d\\ GENROU\\ 2$\r\n",
      "    18,              GENROU e2q 2,                  $e''_q\\ GENROU\\ 2$\r\n",
      "    19,            GENROU delta 3,                 $\\delta\\ GENROU\\ 3$\r\n",
      "    20,            GENROU omega 3,                 $\\omega\\ GENROU\\ 3$\r\n",
      "    21,              GENROU e1q 3,                   $e'_q\\ GENROU\\ 3$\r\n",
      "    22,              GENROU e1d 3,                   $e'_d\\ GENROU\\ 3$\r\n",
      "    23,              GENROU e2d 3,                  $e''_d\\ GENROU\\ 3$\r\n",
      "    24,              GENROU e2q 3,                  $e''_q\\ GENROU\\ 3$\r\n",
      "    25,             TGOV1 LAG_x 0,                $x'_{LAG}\\ TGOV1\\ 0$\r\n",
      "    26,              TGOV1 LL_x 0,                 $x'_{LL}\\ TGOV1\\ 0$\r\n",
      "    27,             TGOV1 LAG_x 1,                $x'_{LAG}\\ TGOV1\\ 1$\r\n",
      "    28,              TGOV1 LL_x 1,                 $x'_{LL}\\ TGOV1\\ 1$\r\n",
      "    29,             TGOV1 LAG_x 2,                $x'_{LAG}\\ TGOV1\\ 2$\r\n",
      "    30,              TGOV1 LL_x 2,                 $x'_{LL}\\ TGOV1\\ 2$\r\n",
      "    31,             TGOV1 LAG_x 3,                $x'_{LAG}\\ TGOV1\\ 3$\r\n",
      "    32,              TGOV1 LL_x 3,                 $x'_{LL}\\ TGOV1\\ 3$\r\n",
      "    33,                EXDC2 vp 0,                     $V_p\\ EXDC2\\ 0$\r\n",
      "    34,              EXDC2 LS_x 0,                 $x'_{LS}\\ EXDC2\\ 0$\r\n",
      "    35,              EXDC2 LL_x 0,                 $x'_{LL}\\ EXDC2\\ 0$\r\n",
      "    36,              EXDC2 LA_x 0,                 $x'_{LA}\\ EXDC2\\ 0$\r\n",
      "    37,               EXDC2 W_x 0,                  $x'_{W}\\ EXDC2\\ 0$\r\n",
      "    38,                EXDC2 vp 1,                     $V_p\\ EXDC2\\ 1$\r\n",
      "    39,              EXDC2 LS_x 1,                 $x'_{LS}\\ EXDC2\\ 1$\r\n",
      "    40,              EXDC2 LL_x 1,                 $x'_{LL}\\ EXDC2\\ 1$\r\n",
      "    41,              EXDC2 LA_x 1,                 $x'_{LA}\\ EXDC2\\ 1$\r\n",
      "    42,               EXDC2 W_x 1,                  $x'_{W}\\ EXDC2\\ 1$\r\n",
      "    43,                EXDC2 vp 2,                     $V_p\\ EXDC2\\ 2$\r\n",
      "    44,              EXDC2 LS_x 2,                 $x'_{LS}\\ EXDC2\\ 2$\r\n",
      "    45,              EXDC2 LL_x 2,                 $x'_{LL}\\ EXDC2\\ 2$\r\n",
      "    46,              EXDC2 LA_x 2,                 $x'_{LA}\\ EXDC2\\ 2$\r\n",
      "    47,               EXDC2 W_x 2,                  $x'_{W}\\ EXDC2\\ 2$\r\n",
      "    48,                EXDC2 vp 3,                     $V_p\\ EXDC2\\ 3$\r\n",
      "    49,              EXDC2 LS_x 3,                 $x'_{LS}\\ EXDC2\\ 3$\r\n",
      "    50,              EXDC2 LL_x 3,                 $x'_{LL}\\ EXDC2\\ 3$\r\n",
      "    51,              EXDC2 LA_x 3,                 $x'_{LA}\\ EXDC2\\ 3$\r\n",
      "    52,               EXDC2 W_x 3,                  $x'_{W}\\ EXDC2\\ 3$\r\n",
      "    53,                   Bus a 0,                    $\\theta\\ Bus\\ 0$\r\n",
      "    54,                   Bus a 1,                    $\\theta\\ Bus\\ 1$\r\n",
      "    55,                   Bus a 2,                    $\\theta\\ Bus\\ 2$\r\n",
      "    56,                   Bus a 3,                    $\\theta\\ Bus\\ 3$\r\n",
      "    57,                   Bus a 4,                    $\\theta\\ Bus\\ 4$\r\n",
      "    58,                   Bus a 5,                    $\\theta\\ Bus\\ 5$\r\n",
      "    59,                   Bus a 6,                    $\\theta\\ Bus\\ 6$\r\n",
      "    60,                   Bus a 7,                    $\\theta\\ Bus\\ 7$\r\n",
      "    61,                   Bus a 8,                    $\\theta\\ Bus\\ 8$\r\n",
      "    62,                   Bus a 9,                    $\\theta\\ Bus\\ 9$\r\n",
      "    63,                   Bus v 0,                         $V\\ Bus\\ 0$\r\n",
      "    64,                   Bus v 1,                         $V\\ Bus\\ 1$\r\n",
      "    65,                   Bus v 2,                         $V\\ Bus\\ 2$\r\n",
      "    66,                   Bus v 3,                         $V\\ Bus\\ 3$\r\n",
      "    67,                   Bus v 4,                         $V\\ Bus\\ 4$\r\n",
      "    68,                   Bus v 5,                         $V\\ Bus\\ 5$\r\n",
      "    69,                   Bus v 6,                         $V\\ Bus\\ 6$\r\n",
      "    70,                   Bus v 7,                         $V\\ Bus\\ 7$\r\n",
      "    71,                   Bus v 8,                         $V\\ Bus\\ 8$\r\n",
      "    72,                   Bus v 9,                         $V\\ Bus\\ 9$\r\n",
      "    73,                    PV p 0,                          $p\\ PV\\ 0$\r\n",
      "    74,                    PV q 0,                          $q\\ PV\\ 0$\r\n",
      "    75,                    PV p 1,                          $p\\ PV\\ 1$\r\n",
      "    76,                    PV q 1,                          $q\\ PV\\ 1$\r\n",
      "    77,                    PV p 2,                          $p\\ PV\\ 2$\r\n",
      "    78,                    PV q 2,                          $q\\ PV\\ 2$\r\n",
      "    79,                 Slack p 0,                       $p\\ Slack\\ 0$\r\n",
      "    80,                 Slack q 0,                       $q\\ Slack\\ 0$\r\n",
      "    81,               GENROU Id 0,                    $I_d\\ GENROU\\ 0$\r\n",
      "    82,               GENROU Iq 0,                    $I_q\\ GENROU\\ 0$\r\n",
      "    83,               GENROU vd 0,                    $V_d\\ GENROU\\ 0$\r\n",
      "    84,               GENROU vq 0,                    $V_q\\ GENROU\\ 0$\r\n",
      "    85,               GENROU tm 0,                 $\\tau_m\\ GENROU\\ 0$\r\n",
      "    86,               GENROU te 0,                 $\\tau_e\\ GENROU\\ 0$\r\n",
      "    87,               GENROU vf 0,                    $v_f\\ GENROU\\ 0$\r\n",
      "    88,             GENROU psid 0,                 $\\psi_d\\ GENROU\\ 0$\r\n",
      "    89,             GENROU psiq 0,                 $\\psi_q\\ GENROU\\ 0$\r\n",
      "    90,            GENROU psiaq 0,              $\\psi_{aq}\\ GENROU\\ 0$\r\n",
      "    91,            GENROU psiad 0,              $\\psi_{ad}\\ GENROU\\ 0$\r\n",
      "    92,             GENROU psia 0,               $\\psi_{a}\\ GENROU\\ 0$\r\n",
      "    93,               GENROU Se 0,        $S_e(|\\psi_{a}|)\\ GENROU\\ 0$\r\n",
      "    94,               GENROU Id 1,                    $I_d\\ GENROU\\ 1$\r\n",
      "    95,               GENROU Iq 1,                    $I_q\\ GENROU\\ 1$\r\n",
      "    96,               GENROU vd 1,                    $V_d\\ GENROU\\ 1$\r\n",
      "    97,               GENROU vq 1,                    $V_q\\ GENROU\\ 1$\r\n",
      "    98,               GENROU tm 1,                 $\\tau_m\\ GENROU\\ 1$\r\n",
      "    99,               GENROU te 1,                 $\\tau_e\\ GENROU\\ 1$\r\n",
      "   100,               GENROU vf 1,                    $v_f\\ GENROU\\ 1$\r\n",
      "   101,             GENROU psid 1,                 $\\psi_d\\ GENROU\\ 1$\r\n",
      "   102,             GENROU psiq 1,                 $\\psi_q\\ GENROU\\ 1$\r\n",
      "   103,            GENROU psiaq 1,              $\\psi_{aq}\\ GENROU\\ 1$\r\n",
      "   104,            GENROU psiad 1,              $\\psi_{ad}\\ GENROU\\ 1$\r\n",
      "   105,             GENROU psia 1,               $\\psi_{a}\\ GENROU\\ 1$\r\n",
      "   106,               GENROU Se 1,        $S_e(|\\psi_{a}|)\\ GENROU\\ 1$\r\n",
      "   107,               GENROU Id 2,                    $I_d\\ GENROU\\ 2$\r\n",
      "   108,               GENROU Iq 2,                    $I_q\\ GENROU\\ 2$\r\n",
      "   109,               GENROU vd 2,                    $V_d\\ GENROU\\ 2$\r\n",
      "   110,               GENROU vq 2,                    $V_q\\ GENROU\\ 2$\r\n",
      "   111,               GENROU tm 2,                 $\\tau_m\\ GENROU\\ 2$\r\n",
      "   112,               GENROU te 2,                 $\\tau_e\\ GENROU\\ 2$\r\n",
      "   113,               GENROU vf 2,                    $v_f\\ GENROU\\ 2$\r\n",
      "   114,             GENROU psid 2,                 $\\psi_d\\ GENROU\\ 2$\r\n",
      "   115,             GENROU psiq 2,                 $\\psi_q\\ GENROU\\ 2$\r\n",
      "   116,            GENROU psiaq 2,              $\\psi_{aq}\\ GENROU\\ 2$\r\n",
      "   117,            GENROU psiad 2,              $\\psi_{ad}\\ GENROU\\ 2$\r\n",
      "   118,             GENROU psia 2,               $\\psi_{a}\\ GENROU\\ 2$\r\n",
      "   119,               GENROU Se 2,        $S_e(|\\psi_{a}|)\\ GENROU\\ 2$\r\n",
      "   120,               GENROU Id 3,                    $I_d\\ GENROU\\ 3$\r\n",
      "   121,               GENROU Iq 3,                    $I_q\\ GENROU\\ 3$\r\n",
      "   122,               GENROU vd 3,                    $V_d\\ GENROU\\ 3$\r\n",
      "   123,               GENROU vq 3,                    $V_q\\ GENROU\\ 3$\r\n",
      "   124,               GENROU tm 3,                 $\\tau_m\\ GENROU\\ 3$\r\n",
      "   125,               GENROU te 3,                 $\\tau_e\\ GENROU\\ 3$\r\n",
      "   126,               GENROU vf 3,                    $v_f\\ GENROU\\ 3$\r\n",
      "   127,             GENROU psid 3,                 $\\psi_d\\ GENROU\\ 3$\r\n",
      "   128,             GENROU psiq 3,                 $\\psi_q\\ GENROU\\ 3$\r\n",
      "   129,            GENROU psiaq 3,              $\\psi_{aq}\\ GENROU\\ 3$\r\n",
      "   130,            GENROU psiad 3,              $\\psi_{ad}\\ GENROU\\ 3$\r\n",
      "   131,             GENROU psia 3,               $\\psi_{a}\\ GENROU\\ 3$\r\n",
      "   132,               GENROU Se 3,        $S_e(|\\psi_{a}|)\\ GENROU\\ 3$\r\n",
      "   133,              TGOV1 pout 0,                 $P_{out}\\ TGOV1\\ 0$\r\n",
      "   134,              TGOV1 wref 0,            $\\omega_{ref}\\ TGOV1\\ 0$\r\n",
      "   135,              TGOV1 pref 0,                 $P_{ref}\\ TGOV1\\ 0$\r\n",
      "   136,                TGOV1 wd 0,            $\\omega_{dev}\\ TGOV1\\ 0$\r\n",
      "   137,                TGOV1 pd 0,                     $P_d\\ TGOV1\\ 0$\r\n",
      "   138,              TGOV1 LL_y 0,                  $y_{LL}\\ TGOV1\\ 0$\r\n",
      "   139,              TGOV1 pout 1,                 $P_{out}\\ TGOV1\\ 1$\r\n",
      "   140,              TGOV1 wref 1,            $\\omega_{ref}\\ TGOV1\\ 1$\r\n",
      "   141,              TGOV1 pref 1,                 $P_{ref}\\ TGOV1\\ 1$\r\n",
      "   142,                TGOV1 wd 1,            $\\omega_{dev}\\ TGOV1\\ 1$\r\n",
      "   143,                TGOV1 pd 1,                     $P_d\\ TGOV1\\ 1$\r\n",
      "   144,              TGOV1 LL_y 1,                  $y_{LL}\\ TGOV1\\ 1$\r\n",
      "   145,              TGOV1 pout 2,                 $P_{out}\\ TGOV1\\ 2$\r\n",
      "   146,              TGOV1 wref 2,            $\\omega_{ref}\\ TGOV1\\ 2$\r\n",
      "   147,              TGOV1 pref 2,                 $P_{ref}\\ TGOV1\\ 2$\r\n",
      "   148,                TGOV1 wd 2,            $\\omega_{dev}\\ TGOV1\\ 2$\r\n",
      "   149,                TGOV1 pd 2,                     $P_d\\ TGOV1\\ 2$\r\n",
      "   150,              TGOV1 LL_y 2,                  $y_{LL}\\ TGOV1\\ 2$\r\n",
      "   151,              TGOV1 pout 3,                 $P_{out}\\ TGOV1\\ 3$\r\n",
      "   152,              TGOV1 wref 3,            $\\omega_{ref}\\ TGOV1\\ 3$\r\n",
      "   153,              TGOV1 pref 3,                 $P_{ref}\\ TGOV1\\ 3$\r\n",
      "   154,                TGOV1 wd 3,            $\\omega_{dev}\\ TGOV1\\ 3$\r\n",
      "   155,                TGOV1 pd 3,                     $P_d\\ TGOV1\\ 3$\r\n",
      "   156,              TGOV1 LL_y 3,                  $y_{LL}\\ TGOV1\\ 3$\r\n",
      "   157,              EXDC2 vout 0,                 $v_{out}\\ EXDC2\\ 0$\r\n",
      "   158,                EXDC2 Se 0,                     $S_e\\ EXDC2\\ 0$\r\n",
      "   159,              EXDC2 vref 0,                 $V_{ref}\\ EXDC2\\ 0$\r\n",
      "   160,                EXDC2 vi 0,                     $V_i\\ EXDC2\\ 0$\r\n",
      "   161,              EXDC2 LL_y 0,                  $y_{LL}\\ EXDC2\\ 0$\r\n",
      "   162,               EXDC2 W_y 0,                   $y_{W}\\ EXDC2\\ 0$\r\n",
      "   163,              EXDC2 vout 1,                 $v_{out}\\ EXDC2\\ 1$\r\n",
      "   164,                EXDC2 Se 1,                     $S_e\\ EXDC2\\ 1$\r\n",
      "   165,              EXDC2 vref 1,                 $V_{ref}\\ EXDC2\\ 1$\r\n",
      "   166,                EXDC2 vi 1,                     $V_i\\ EXDC2\\ 1$\r\n",
      "   167,              EXDC2 LL_y 1,                  $y_{LL}\\ EXDC2\\ 1$\r\n",
      "   168,               EXDC2 W_y 1,                   $y_{W}\\ EXDC2\\ 1$\r\n",
      "   169,              EXDC2 vout 2,                 $v_{out}\\ EXDC2\\ 2$\r\n",
      "   170,                EXDC2 Se 2,                     $S_e\\ EXDC2\\ 2$\r\n",
      "   171,              EXDC2 vref 2,                 $V_{ref}\\ EXDC2\\ 2$\r\n",
      "   172,                EXDC2 vi 2,                     $V_i\\ EXDC2\\ 2$\r\n",
      "   173,              EXDC2 LL_y 2,                  $y_{LL}\\ EXDC2\\ 2$\r\n",
      "   174,               EXDC2 W_y 2,                   $y_{W}\\ EXDC2\\ 2$\r\n",
      "   175,              EXDC2 vout 3,                 $v_{out}\\ EXDC2\\ 3$\r\n",
      "   176,                EXDC2 Se 3,                     $S_e\\ EXDC2\\ 3$\r\n",
      "   177,              EXDC2 vref 3,                 $V_{ref}\\ EXDC2\\ 3$\r\n",
      "   178,                EXDC2 vi 3,                     $V_i\\ EXDC2\\ 3$\r\n",
      "   179,              EXDC2 LL_y 3,                  $y_{LL}\\ EXDC2\\ 3$\r\n",
      "   180,               EXDC2 W_y 3,                   $y_{W}\\ EXDC2\\ 3$\r\n"
     ]
    }
   ],
   "source": [
    "!cat kundur_full_out.lst"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Windows users may not have `cat` by default in the search path. Do instead\n",
    "\n",
    "```\n",
    "with open('kundur_full_out.lst', 'r') as f:\n",
    "    print(f.read())\n",
    "```"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Plot and save to file"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "ANDES 0.6.9 (Git commit id c8cba088, Python 3.7.6 on Darwin)\n",
      "Session: hcui7, 02/23/2020 02:16:40 PM\n",
      "This program comes with ABSOLUTELY NO WARRANTY.\n",
      "\n",
      "Figure saved to <kundur_full_out_1.png>\n",
      "Figure(640x480)\n"
     ]
    }
   ],
   "source": [
    "!andes plot kundur_full_out.lst 0 2 8 --save"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Display image"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAoAAAAHgCAYAAAA10dzkAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAPYQAAD2EBqD+naQAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAgAElEQVR4nOzdeXxbV5028EeLV3mR5S1x7MSWm7RZ2rTOAi1pKVQutFC22ilLedlaCRg6zAC18cwAAwMYi3mZGWYYkDJlZ0piQYEWaGvx0n2zrSZdsktO4ix2bEtX8iLJsnXfPxypli3Jkn1l2dHz/Xz8Aenec+5P1EWPz7nnXJkoiiKIiIiIKGPI010AERERES0vBkAiIiKiDMMASERERJRhGACJiIiIMgwDIBEREVGGYQAkIiIiyjAMgEREREQZhgGQiIiIKMMwABIRERFlGAZAIiIiogzDAEhERESUYRgAiYiIiDIMAyARERFRhmEAJCIiIsowDIBEREREGYYBkIiIiCjDMAASERERZRgGQCIiIqIMwwBIRERElGEYAImIiIgyDAMgERERUYZhACQiIiLKMAyARERERBmGAZCIiIgowzAAEhEREWUYBkAiIiKiDMMASERERJRhGACJiIiIMgwDIBEREVGGYQAkIiIiyjAMgEREREQZhgGQiIiIKMMwABIRERFlGAZAIiIiogzDAEhERESUYRgAiYiIiDIMAyARERFRhmEAJCIiIsowDIBEREREGYYBkIiIiCjDMAASERERZRgGQCIiIqIMwwBIRERElGEYAImIiIgyDAMgERERUYZhACQiIiLKMAyARERERBmGAZCIiIgowzAAEhEREWUYBkAiIiKiDMMASERERJRhGACJiIiIMgwDIBEREVGGYQAkIiIiyjAMgEREREQZhgGQiIiIKMMo013AahYMBnH+/HkUFhZCJpOluxwiIiJKgCiKGB0dRVVVFeTyzBwLYwBcgvPnz6OmpibdZRAREdEi9Pf3o7q6Ot1lpAUD4BIUFhYCmPkFKioqSnM1RERElAiPx4Oamprw93gmYgBcgtC0b1FREQMgERHRKpPJt29l5sQ3ERERUQZjACQiIiLKMAyARERERBmGAZCIiIgow3ARCBERLZooipiensbU1FS6SyGCUqmEQqHI6MUdiWIAJCKipImiCEEQMDQ0hOnp6XSXQxSmUChQUVGB4uJiBsE4GACJiChpAwMDEAQhvA2WUqnkly2llSiKmJqagsfjwYULF+D1erF27dp0l7ViMQASEVFSpqen4Xa7UV5ejrKysnSXQxShsLAQOTk5GB4eRkVFBRQKRbpLWpG4CISIiJISCAQgiiJUKlW6SyGKSqVSQRRFBAKBdJeyYjEAEhHRonDKl1Yq/m4ujAGQiIiIKMMwABIRERFlGC4CoRXhoseHoAisKc5NdylERJKyWCzo6uqCWq0GAJSWlkKv16OnpwcAoNPpwufabDaYTCaYzWY0NTVh165d4WN2ux0HDhxAR0cH9Hp91DYtLS3o6OgIHxMEAffeey+sVit0Oh327dsXriPWtex2O9RqdUQ/c5nNZgiCEO7Lbrejra0t/Hpu/zqdDs3NzeG6HQ4HWltbYbVaodfr57WlZSDSorndbhGA6Ha7013KqiVMTIrtfzoibvzHP4l3/OfT6S6HiBLg9XrFw4cPi16vN92lrGgul0tsamoSOzo65h1raWkR1Wp11Ha9vb0iANHlcs071tnZKXZ2dkZto9fro/bpcrmi1jD7WnPp9XpRp9PNe99ut4s6nU7s7e2ddw2tVjuv5nifJV5dc3V0dIgmk0k0mUwJtVnod5Tf36LIKWBKi2BQxAPP9OGt3/0rfvbcKexYX4LXz3vgneSGskS0+jkcDtTV1cFgMKClpWXe8cbGRmi12qhtrVYrGhoaoo6IabXaqO16enrQ0dEBQRBgtVrnHZs9YhjtWnPt2LFjXj+CIKCxsREmk2leG7VaDYPBgHvvvXde/1qtNupnsVqtaGpqilrXbEajEWq1Gnq9Hnq9Hg0NDTAYDAu2o/gYACktHj88gH955DBu27YGT95/M/7p3ZsxHRTx2nl3uksjIlqyxsZGtLW1RUzvzrZz586Yx7q6urBz587wa0EQIAhC+HWs4KhWq9HU1ASTyRTxvsPhiDm92tXVFbWO3t7eeSGvubkZBoMh5vW1Wu280Bir/1Bdsfqarb29HXv37g2/1ul0MJvNC7aj+BgAKS2ODYyhVJWN9g9cg4qiXFxZWYjcLDkO9QsLNyYiSqG5gWs2h8MxL+TMZTQa4XQ6o478hYRGzKKxWq1obGwMvzabzeEAF2tkMMRgMMBiscSsf6Frhd6zWq3o7OwMv2exWNDT0xP3MzmdznnXjdZ/MhwOR8S9hnP7psXjIhBKi9POcWwozQ+/VirkuHpdMV5mACRa1byT07APjaW7DABAfXkB8rKTfwpEaDFFNFarFSaTCb29vTHbm0ymiBGrWKKNftlsNgBAd3c3HA4H7HY7HA5H3ODlcDjCI4Y6nQ5qtTruZ5h7LYfDAYvFAofDgZGREZSWlqK3tzcidLW3t8ccyQsJLR6Z23+0drFCXbTPFo1arU445FJ0DICUFqdHJrChNPIpAtur1Xj09YE0VUREUrAPjeHd//lMussAADxy3x5sW1ecVBtBEOJOS+r1enR1dcVt73A4sGPHjqSuG9LT0wOtVhtegWuz2bB///64bUIraWfXaDKZ0NLSApvNFjO4he7/C7UNrRiur6+PCGeCIMBms+Guu+5asI7Z14p3/19PT8+CgTIejUYDp9O56PbEAEhpcnpkAnuuiHyG6PYaNf7nmT4Mj/lRVpCTpsqIaCnqywvwyH170l0GgJlakjV3YYLNZoPT6YwIK/GmNEOhJFqItFgs6O7uhtFoRENDA9ra2uYtgujq6oq4906r1UZcL5H75gwGA4xGI2w2W9wFIHPvzwtNSzc2Nka0CX2maItFQmw2G2w2W8TWMd3d3TFDXrxgmgiGv6VjAKRlN+6fwvCYH7Vl+RHvX1sz81fioX4Bt2yuTEdpRLREedmKpEfdVrK5o2uzp1ujCYWzaNOTTU1N0Ol0MBqN6OjoiBqArFZrRIhSq9Xh8xwOR0IBUKvVoqGhASaTKe5IZE9Pz7z7EKNNuSayUKO1tRU6nW7enoatra0Lto0n1rUXGqmlhXERCC270yMTAID1msgp4OqSPJSqsrkQhIjSpru7O+L1/v37I6Ywo22BMldTU9O8fkJCCxeihb/QgodYI2Mmk2nesVhByGAw4MCBAzFrjHWt3t7eqP01NDSE7+mbK3T/4OxFI8BMeNVoNPPON5vNMUcl5wpNIUcLpksZQSQGQEqDM85xAEBtaeQIoEwmw/YaNQ6e5VYwRJQeoX30bDYbmpubodVqYTQa4XA4YDAYIp7MEcu+fftgNpvnBSZBELB///6YAdJqtUKtVsecPo7GbDZHHZHU6/UQBCHmaKXFYol6f57T6Qy/Z7VawyOZ+/btm7e9TKifrq6ueQtAgJkQOvd+yVjXjaetrS1ixa/FYkk4QFJsnAKmZXdqZAIFOUpoVNlvvPnoPwBjA7i2pg0PPNMHURQhk8nSVyQRZSSDwYBbbrkFGo0GJpMJO3fuxC233IL29nbo9fqENi5Wq9Xo6+tDa2sr1Go1SktLw4Gns7Mz6vYlra2t4ZBnNBrD74+MjIQDqd1uD78vCAJaW1tx4MABdHd3RzziLaSlpWVe2BQEAe3t7eHRP6PRGLFauK2tDe3t7eFaZm8/09nZCYPBEJ5WDo0+RguGwEwINZvNaG1tRWlpKYCZUbuFRlDnamlpgdFoDO/9Z7fbY16TEicTRVFMdxGrlcfjQXFxMdxuN4qKitJdzqrR9ttXcahfwJ8+f+PMG8Eg8I0SAMDz73saH/p1P/76pZtRV6aK0wsRpYvP50NfXx/q6uqQm8vnd9PKs9DvKL+/OQVMaXDGOR65AEQuB758BlDm4lrXzHQB7wMkIiJKHQZAWnanhifeWABy7FHA9gsgpwi48jbkHfkN6spUOMgASERElDIMgLSs/FPTuOD2ziwACXiBP90PHH0EkMmAnZ8ErnwnrlvHAEhERJRKDIC0rM66vAiKwPrSfOD5/wJGLwC3fmvmYN1NwC1fxdXry3D4vAf+qen0FktERHSZYgCkZXXm0h6AdTmjwNPfA978aaDsijdOGB3ErcIBTE1P4eiF0TRVSUREdHljAKRldXpkHNlKOSovPAFM+YCb7o88wd2Pdd3fxo3Kwzh0ltPAREREqcAASMvq1MgEakryIN++F/j0M0DunEdGrdsBaOpxd97zsF8cS0+RRERElzkGQFpWZ5wTqC1VAdkqoHLr/BNkMuCau3Dj1AsYGB5Z/gKJiIgyAAMgLatTI+MzC0D+8LfAkYejn3TNXuSKXtQMPbm8xREREWUIBkBaNtNBEWedXtSqs4CDvwJGB6KfqKlD98bP45WxYgSDfFANERGR1BgAadkMeHyYnA7iyjwBCE4BZRtjnutp+Bu8NHUFLo76l7FCIiLpWSwWGAwGtLa2orW1FUajEYIgwGq1znsusM1mg8FggEwmQ3NzM4xGY/jHYDCgpKQk/EzcaG1aW1sjjgmCgObmZpSUlKC5uRmCICx4rVCt8ZjN5vDzeUPP+53d99z+GxsbI+p2OBzhuqK1pWUg0qK53W4RgOh2u9Ndyqrw7IkhcUPrI+L5lx4Sxa8ViaLQH/Ncx8mj4t+2fVl80TGyjBUSUSK8Xq94+PBh0ev1pruUFc3lcolNTU1iR0fHvGMtLS2iWq2O2q63t1cEILpcrnnHOjs7xc7Ozqht9Hp91D5dLlfUGmZfay69Xi/qdLp579vtdlGn04m9vb3zrqHVaufVHO+zxKsr2nVNJlPUmqJZ6HeU39+iqExr+pzF4XDAarWis7MTXV1dCbUxGo1Qq9UAZv7KaWlpSeo4ALS2tqK+vh4AoNFo0NTUtJSPQXGcdk5ALgPK/f1AVj5QWBXz3Krxw/iP7P/GIxc+CNRplrFKIqKlczgc2LFjBzo7O6HT6eYdb2xsnDf6F2K1WtHQ0BD+/ppNq9VGbdPT04OOjg6YzWZYrdaIa/b09ECv18e91lw7duyYN9IoCAIaGxvR1dU1rw61Wg2DwYB7770XnZ2dEf1rtdqon8VqtSb0nWuz2dDT0wNBEOB0Ohc8nxKzIqaAbTYbrFZrUv9wQ+FOr9dDr9ejoaEBBoMh4eOCIGDHjh1oa2uDXq/Hzp070dzcLPlnozecGhlHlToPyrobgFu/Cchj//rlrN0CAPCee325yiMikkxjYyPa2tqihj8A2LlzZ8xjXV1d2LlzZ/i1IAgRU6SxQqBarUZTUxNMJlPE+w6HI2oAC10rWh29vb3zgmFzczMMBkPM62u12nmhNlb/obpi9TVbQ0MD9Hp9QudS4lZEAFzMP9z29nbs3bs3/Fqn00X8tZLI8bvuuiv8L0VDQ0PCI4+0OGdGJrChNH9mr79dn4p/skaLAJRQDB9dnuKIiC6ZG7hmC81WxWM0GuF0OqPOOoWERsyisVqtaGxsDL82m80R31WxwhwAGAwGWCyWhO+pm3ut0HuhGbkQi8WCnp6euJ/J6XTOu260/mllWDFTwMlwOBwQBCHmkLJWq417XKfTwWg0ore3Fw6HAw6HAzqdLuZfKSSN0yMTuHa9GrD9AthwA1BaH/tkRRaGcmqg8pxYvgKJSBqjA/NX+eepgZJaIOADhqL8YVd17cx/Dp8AJscjj6nXA/kaYHwYcJ+NPJZTOPP/JcFpYODVN94vXDPzswhmszlm0LFarTCZTOjt7Y3Z3mQyRQxAxBJt0MNmswEAuru74XA4YLfb4XA44gYvh8MRHjHU6XRQq9VxP8PcazkcDlgsFjgcDoyMjKC0tBS9vb0R36Ht7e0Lfkfa7faINqH+o7WL9R1Ny2fVBsBo1Go1BEFI+Hho+Fmr1cJgMKC5uTnuL7jf74ff/8aqVI/Hs4RPkVlEUcTpkXHcuU0N/OFzwPtN8QMggPOlb8GZAa4CJlp1en4CPPmdyPeu3gvcuQ/wnAPMb53f5p/dM//5u88AZ7sjj73fDGy/C3j9IeBPX4o8Vv924KMPzYTG2f2+9cvA29qSLl0QhLizUXq9Pu5sUeg7ZseOHUlfG5i5X0+r1aKjowPATIjav39/3DZWqzXiHj+9Xg+TyYSWlhbYbLaY32uh+/9CbQVBwL333ov6+vqIcCYIAmw2G+66664F65h9rXj3//X09HDQJc1WZQCMRaPRwOl0xvyrInQ8FADVanX4HoeOjg7U1dXB5XLF7L+9vR1f//rXpS88A7i9AYxPTmNT1sWZN0qvWLDNqR1t+FbnIXw0MI3cLEWKKyQiyez8BHDlbZHv5V36/+WidYA+zibv7/th9BFAANj6fqB6V+SxnMKZ/8xWRfa7yNG/uQsTbDYbnE5nRFiJN6UZuo89Woi0WCzo7u6G0WhEQ0MD2tra5i2C6Orqirj3TqvVRlwvkfvmDAYDjEZjePFErAUgc+/PC01LNzY2RrQJfaZoi0VCbDYbbDZbOLgCM6OYsUJevGBKy2NF3AMolYUWkISOazQzq0pn32QbGh2Md29HW1sb3G53+Ke/v1+CqjODc3wSAFA5eel/swVG/wBgvSYfJfDg7HDsUE5EK1Dhmpkp3dk/JbUzx7Jy5x8LTf8CM/uDzj2Wf2knAFXZ/GOh/y+RKyLfX2QAnMtqtUZ8V8yebo0mFM6i3YPX1NSEtraZUcmOjo6oK2Dn3jOnVqvDQSl0y9JCtFotGhoa5i0Gmaunpwe7dkUG6mj9J3J/fmtr67xbqWw226JHQin1VmUAjPXLGBq6T/T43MCoVqvj/suVk5ODoqKiiB9KjGsiAAAo8Z0B8suAvJIF29TJLuDl3E/Dc/zZVJdHRARgZtRqtv3790fMKplMprgjYcBM0JvbT0hokCHa6Ffo/vZYI2Mmk2nesVhT1gaDAQcOHIhZY6xr9fb2Ru2voaEhfE/fXKH7B2cvGgFmvlNDAy6zmc3mmKOStHxWbQCMFdZ0Ot2Cx9VqNbRa7bzjgiDE/cuOFk+YmBkBzCmuBDa/O6E2pdWbMCkqMXn+cCpLIyIKC80E2Ww2NDc3Q6vVwmg0wuFwwGAwzBsxi2bfvn0wm83zApMgCNi/f3/MAGm1WsPfT3NZLJaobcxmc9TvLb1eH/c7zWKxRL0/b/ZtVKHt2UKfKdqIosViQVdX17wFIMBMCJ17v2Ss6yaCewBKa0XdAxjrH25ohdLsFU1tbW0RN75aLJaIvygWOt7R0RFx/4PFYoFOp1vwLztanNAIYM719wDKxO7nkyuz0K9YB+XIsVSWRkQUZjAYcMstt0Cj0cBkMmHnzp245ZZb0N7eDr1en9DGxWq1Gn19fWhtbYVarUZpaWk48HR2dka91ai1tTUc8oxGY/j9kZGRcCC12+3h9wVBQGtrKw4cOIDu7m7s27dvXqhqaWmZ950mCALa29vDo39Go3Hed2t7e3u4ltnbz3R2dsJgMISndUOjj7GmmvV6ffgxcaWlpQCwqO9Zh8MBk8kU/t8h9AAHjiIujUwURTHdRcz9h9vS0hLxD9dsNqOjoyPilx+IfNKH3W6PuPk0keNmszn8183IyMi84wvxeDwoLi6G2+3mdPAC/udpB77XdQyHv7AVKKqOuwn0bC8Y3wvN9Ag2tT2T4gqJKFE+nw99fX2oq6tDbm5uusshmmeh31F+f6+QALha8Rcocd997CietB3GI/5PAHt/AWx5T0Lt/mL6Eq4ZeAjlX7MvfDIRLQsGQFrpGAAXtirvAaTVxzURwJbswZkXZRsTbnfmyk/i5qn/BP9OISIikg4DIC0LYWISm5SDAGSAJvFH/lWVazAeEDFyaRsZIiIiWjoGQFoWrvEAanFhZkNXZU7C7daX5OHX2f+C8ed/nMLqiIiIMsuKWgVMly/XxCQqZCMJPQFktppSFYBxBM9F33+KiIiIkscRQFoWrolJWLd8G/jgr5JqV5CjxGnFeuQ4j6eoMiIioszDAEgpJ4oiXBMBlORnAVl5SbcfztNCPW4HuBCEaEXh4ixaqfi7uTAGQEo5b2AaxVNOfOCFZuBsT/Lt1RuRPz0KjF1MQXVElKysrCzIZDKMj4+nuxSiqMbHxyGTyZCVlZXuUlYs3gNIKeeaCKBC5kLx6AlAlvzfHL6qN+Fvh76G7+dm5l5NRCuNQqFAcXExhoaG4Pf7UVRUBKVSCZlMlu7SKIOJooipqSl4PB54PB6o1WooFIk9eSoTMQBSyrnGJ1Euc8+8KKhIun1F5Vp8b+xK/KssB9kS10ZEi7NmzRrk5eXh4sWL8Hg86S6HKEyhUGDt2rUoLi5OdykrGgMgpZwwEUBZKACqypNuX6PJx/+RPwZXtweV139I4uqIaDFkMhnUajWKi4sxPT2NqampdJdEBKVSCYVCwdHoBDAAUsq5JiZRBjfEnGLIktgDMKSmJB/vlHdDdmQIYAAkWlFkMhmUSiWUSn6dEK0mXARCKSdMTKJL3A18wLSo9muLc3EGlVAIp6QtjIiIKEMxAFLKuSYCcOdvgOzK2xbVXqmQw5WzDqrxfokrIyIiykwMgJRyrolJfFD5BHCia9F9eAvWI3d6FPC6pCuMiIgoQzEAUsq5xiexd+ph4MTji+7DW7oVD+e/HwhOS1gZERFRZmIApJRzTQSgDgqAKvktYEJy12zCN6fuBlRlElZGRESUmRgAKeU8414UBN1AQfJbwIRUq/OgGT2OyYHDElZGRESUmRgAKeXEiRHIIS5pBLC6JA/tWfvgf/I/JKyMiIgoMzEAUsq5vZM4XnkbUFq/6D6qS/LRL1YgONInYWVERESZiQGQUmpqOog+XyFe3mkEyq9cdD9rinNxRqyE0nNawuqIiIgyEwMgpZTgDaAEHqwVB5fUT7ZSDiG3Gvm+QWDKL1F1REREmYkBkFJKmJhEk+IpvKXrfUvuy1dUi6HsamB8WILKiIiIMhcDIKWUayKAMpkb0/mLXwEcMla5G3+jMQPF6ySojIiIKHMxAFJKucYnUSZzQ1aw+BXAIdUleTjr8gLBoASVERERZS4GQEop18QkyuGGorByyX2tU+fh+942TD/aJkFlREREmYsBkFLKNRFAoWIS8kIpRgDz4RHz4b94UoLKiIiIMhcDIKWUa2ISn8/vAG777pL7qi7JwxmxAnCdWnphREREGYwBkFJKGA+gJD8LkC/9V22tOhenxUrkjJ7hfYBERERLwABIKRUYHcIP3H8DnH5+yX3lKBXw5K2DIjgJjF6QoDoiIqLMxABIKSUbG0R14BQgV0jS32DJTny7/ldA4RpJ+iMiIspEynQXQJc3hXdo5r+olr4PIACUlZbioNsnWaAkIiLKRBwBpJTK9o3M/BcJ9gEEZraCee/FHwG2n0vSHxERUSZiAKSUEUUR+ZMjCCjygGyVJH1Wl+SjPnAMQftfJemPiIgoEzEAUsqM+afw5+kdsO3+nmR9Vpfk4XSwEoFhh2R9EhERZRoGQEoZ13gA/WIlAtpbJesztBegwtUnWZ9ERESZhgGQUsY1MYlmxROoGXpSsj6r1HnoE9dAOekGxkck65eIiCiTMABSyrgmJvFxxWMou/CEZH3mZilwJn8rnqj7AlcCExERLRIDIKWMMBFAmcyN7OJKSfvN0tTg4bz3AXlqSfslIiLKFAyAlDKucR9K4YGySNoAuE6dB83A00Df05L2S0RElClWzEbQDocDVqsVnZ2d6OrqSqiN0WiEWj0zCiQIAlpaWhI+brVaYTKZ0NjYCK1Wi66uLuzatQtNTU0SfSLyekaglAUl2wMwpLokH7ecPAC89AJQd6OkfRMREWWCFTECaLPZYLVaIQgCnE5nQm1C4U6v10Ov16OhoQEGgyHh44IgwGazwWAwwGAwoL6+nuFPYqMTfliz3gaUbpS03/WafBwNVCA4fFLSfomIiDKFTBRFMd1FhFgsFrS3t6O3t3fBc0tKStDX1xce4QMAmUyG0MdZ6LjFYoFOp4s4niyPx4Pi4mK43W4UFRUtup/L1X0PvozhUT8e1L9Z0n6ft4/g8R9/DV/JPQD5P17gYhAiIkoKv79XyAhgshwOBwRBiBrerFbrgsdpeUyNDmFj1kVA4r8xtOUqOMQqyKf9gPuspH0TERFlghVzD2AyHI7oT4FQq9UQBGHB4yEHDhyARqMBAHR3d6OjoyPudf1+P/x+f/i1x+NJtvSM0iA8jo95fw5gUNJ+KwpzcE5Zg7Mlb0b1lE/SvomIiDLBqhwBjEWj0cS9h3D2ca1WC51Oh6amJjQ1NaG+vh7Nzc1x+29vb0dxcXH4p6amRtL6Lzf5AScmsjSATCZpvzKZDDlltfhBzXeB8isl7ZuIiCgTXFYBcKEFJLOPNzQ0QKvVhl/v3bsXFoslYoRwrra2Nrjd7vBPf3//0ou+jBVMOeHLKU1J33VlKvRfdAKj0o4uEhERZYJVGQBnB7fZBEGAVqtd8DgwswhkttD9grGmjwEgJycHRUVFET8UW/G0gMncspT0rS1T4VOD3wZ+9+mU9E9ERHQ5W7UBUK1WRw1rOp1uweOCIKC5uTnieGjkL1Z4pOT4AtPIhQ/T+akJgHXlqktbwZxISf9ERESXsxUVAGNN4TocDhiNxoj32traIlb0WiwW6PX6hI6r1WqYTKaIsGc2m9HU1LSkbWHoDaO+KXxw8is4seubKem/rqwADnEtZO6zQMCbkmucE7x43j6Skr6JiIjSaUWsAnY4HDCZTLBarbDZbGhtbUV9fX04sIWe2jH7SR4tLS0wGo0wm80AALvdDpPJlPDxvXv3RoTKkZERdHZ2pvRzZpJRXwAAUJiXnZL+60pVcATXQgYRcPYBlVsk69vjC+CHT9jxwDN9CEwH8eI/3IKKwlzJ+iciIkq3FbUR9GrDjSRjO3R6BLkP3ITcd34dG25IzRNW3v4NC/4SvAeyD+8HNr1Dkj7/98Uz+NfHj2Ficgofu74W//NMH77+nq24+80bJOmfiIjSj9/fK2wKmC4f3lEXrpSfRb4idX9flJStReuVj0kW/l4+48I/PPQqbt5Ujie+9EpoJaYAACAASURBVDa03b4Zb9Zq8OhrA5L0T0REtFIwAFJK+EeHAQC5RalZBAIAdeUFOOackqy/R18bQFlBNr7bvB1rimemfN+5bS2ed4zANT4p2XWIiIjSjQGQUiIwOrN4Irc4NfsAAjN7Ad489EuI++9ecl+iKOLR1wfQuGUNFPI3Nq5+x9ZKBEURXYe53yAREV0+GAApJabGZwJglip1AVBbpsJUYApi3zNL7uvowChOj0zgndvWRLxfUZiLXRs0+PNrF5Z8DSIiopWCAZBS4kzORrQo7gdU5Sm7Rl25Cn3iGsh9LmAi/lNgFvLoawO4JncQN/1hD2D9ZyDwxjOGb7t6DZ45OQzPpZXNREREqx0DIKXEYLAYvfl7AGVqtoEBgNpSFfrEtTMvlrgh9GOvD+A7qv2QBQPA8z8ATDcCZ3sBAO/ctgaBaRF/OcJpYCIiujwwAFJKVA6/gDvFrpReIzdLAX9R3cyLkZOL7qdveBxHB0bh2XUf0PRjwPAUkFMETMxMY68tzsO1NWr8+VWuBiYiossDAyClxBWup/Ee38Mpv05VuQb/tfbbwBW3LLqPx187h7wsObbfcBugvRmo2Ax8qgvYdOvMCcEgbr96DZ48PoRxv3SrjomIiNKFAZBSIjvghk9ZnPLr1JWp8Ij3aqBwzcInxyDv/h88rPom8hTBWW/KgekA8IM3Abaf4bZta+GfCuKJY0MSVE1ERJReDICUErkBN/zZqd9dvbZMhbzh1yA+8gUgGFy4wRwDA+fQPPZLKCs3A4qsyIOKLCBPA5x4HDWafGytKsLjhzkNTEREqx8DIKVE3vQoprJTPwKoLVMhP+iBrOcBYPh40u1df/wGFAii9D3fiH7CplsBxxNAwIeG9SU4emF0aQUTERGtAAyAlBLd4mYMlTSk/Dp1ZSq8HNwIUaYAzjyfXGOfG/X9v8Hj6rtQWFoV/ZyNtwKBCeD0s6gvV6FveBzTQT4+m4iIVjcGQEqJ9sAH0V97Z8qvU12Sh0l5HkYKr0o6AE5cOAqnWAjxuo/EPqliC1BUDTieQH1FASangzjrmlhi1UREROnFAEiSmwxMY830eRRnTaf8WkqFHOtL83Ey92rgdHIB8FlvLa73fx+7r94W+ySZDPjknwHdP6O+vAAAYB8aW0LFRERE6ccASJIb8zjxZM4XoHU+tSzX21RRiIen3wzc8DkgmGDonHDi4OHDWF9agPWl+fHPVa8H5AqsKcpFfrYC9ovjSy+aiIgojRgASXIT7mEAQFZB6p4DPNuuOg0sg2vg33EPIFck1qj7AXzmtQ/j5voEViqLIvDz90LebYa2XMURQCIiWvUYAElyXs9MAMxepgC4u1YD/1QQfT2PA0f/uHADUUTg5Qfx2HQDbrhy3cLny2SATA4c+xPqywsYAImIaNVjACTJTY7OPEItr6hsWa63eW0hVNkKBA8+CPy/by3c4FwvsgQ7fhe8CdfXJxhSN74DOPUsrtLIYR/iFDAREa1uDIAkuclxNwAgT12+LNdTKuTYUavBc4GNwMXDgNcVv8GhB+FSlMG/7gYU5WbFPzdk061AMIBdwVfgHJ+Ec3xy6YUTERGlCQMgSe5k2dtR7/sFCgvVy3bN3bUl+M3wegAi0N8d99ygKOLX0zfjLZuSeHycRguUbsQV7ucAAA5OAxMR0SrGAEiSG/VNQanMQnZWggsyJLCrVoMj/lIE8sqBM8/FPfflq7+CDt8HcNOmJKeoP7wfubd/CzIZt4IhIqLVjQGQJLfx5E/w71k/WNZrbq9RI1uhwGtr7wRK6qKfJIrAK5149tgFFOUqcU11kiOUpfXILdSgpiSf9wESEdGqxgBIkisaO4n1sovLes3cLAW21xRjn2IvsONj0U96tRP47T0YPPIM9mwsg0IuS+4irtPArz+C69UC7Bc5AkhERKsXAyBJTjnpxoQigf31JLa7ToNuxwjEcy8DTkfkQf8Y0PVVBDbdgV8PVuPGjYtYoJJdABx9BLtz+zkFTEREqxoDIEkuZ9IDX1bxsl93V60GI+N+BDo/Bfz0DmDE/sbBZ/8dmHDiufrPYzooYs8Vi9iiRlUKFKzBlfJ+nHFOwD+V+kfdERERpQIDIEkud8qNQNbyjwDu2FACyOT48w4zkJ0P/OR2YOgY4D4HPPt94Ib78Ns+JbTlKtRoFnj8WyyVW7DO70BQBE6PTEj7AYiIiJYJAyBJ7sf5n8QrZbcv+3ULc7OwpaoIT15QAh//I5CvAX76rpl9Ad//Q5zebMDDh87j4zfULv4iFVtQ5DkOALwPkIiIVi0GQJLck+K1GC3ZmpZr76rVoPuUEyioAD72CLD2WkCZA2y7E//17AWUFeRg786axV/g6ibI3/5PUOcp4RjmSmAiIlqdGABJWgEf3jtuQbU4kJbL767VoN/pxQW3d+aevbstQNlG9Dsn8NDL56C/SYvcpexPWHUdZNvvQn1FIUcAiYho1WIAJGlNDONvg7/A2uCFtFx+V50GAPDTZ09BFMXw+z980o7ivCx85E0bln6Rg/+Lt6scXAlMRESrFgMgSSowNgIAUOaXpOX6ZQU5aH3nVTA95cDXHz6MYFDEBbcXlp6zuOdGLfKyJXg6yTP/hrf6n4J9aDwiZBIREa0WynQXQJcXn2cEWQCyCzVpq+EzN9ejMFeJr/z+Nbi9AahyFMjPUeCj10sw+gcAFVuwbqgPY/4pXBz1o7IoV5p+iYiIlgkDIEnK5xlGIYCcwkXssyehu9+8AUV5WfjC/oOYCor4QuMmFORI9OteuRVF9icBiLBfHGMAJCKiVYcBkCQ1Ji/CoenrUFGYning2d6zvQqFuUr84vnT+NhStn6Zq2ILFH4X1skF2IfHccNiNpUmIiJKIwZAktSAZhfuCdyPv+avjFGxt11ZgbddWSFtp2u3A1s/gHV2Oc65vNL2TUREtAy4CIQk5XNfRBHGpJtuXYnUNUDzTyDX1OGcwABIRESrDwMgSar24Hfx8+zvoDD3Mg6AADA+jOvyBnGeAZCIiFYhBkCSlMwnwIPCpW22vBr85ev4xMC3OAVMRESrEgMgSUrpFzChKEx3GalXsRWl3lMYHh1HYDqY7mqIiIiSsmLm6RwOB6xWKzo7O9HV1ZVQG6PRCLVaDQAQBAEtLS1JHZ+tsbEx4etSbFmTbngVlekuI/Uqt0IRnMQGDGDA7UONJj/dFRERESVsRYwA2mw2WK1WCIIAp9OZUJtQuNPr9dDr9WhoaIDBYEj4+GwWiwVWq1WSz5LpcqY88GcVpbuM1KvcCgDYLDvDhSBERLTqyMQV9Cwri8WC9vZ29Pb2LnhuSUkJ+vr6wiN8ACCTycKP5lroeIggCDhw4AAMBkPSj/XyeDwoLi6G2+1GUVEGhJ4E3PerbrjHvfi5/qZ0l5Jywe834CuDb0XDB76IO3dUp7scIiJKEL+/V8gIYLIcDgcEQYgIdyFWq3XB47MdOHAAe/fuTVmtmcbjF5GXlxnTofK/teHR3Nu5EpiIiFadVRsAo1Gr1RAEYcHjIVarFTqdLuHr+v1+eDyeiB+aZWwIX77weWwRT6a7kmWzTp3LKWAiIlp1VmUAjEWj0cS9h3DucUEQoNVqE+6/vb0dxcXF4Z+ampol1XvZmRjG5sBhFGSlu5Bl0vtT/NT9SQZAIiJadS6rALjQApLZx81mM5qampLqv62tDW63O/zT39+/qDovW14XAECer0lzIcskpwiaqYsYdQ2muxIiIqKkrJhtYJIRa9QuNKK30HGbzYadO3cmfd2cnBzk5OQk3S5jXAqAWQUZEgBL6wEAOe5TEEURMpkszQURERElZtUGQLVaDYfDMS/she7pi3fcarWGt54BALvdDmBm6xitVpv0yCDNCE44IQeQXVCS7lKWh2bmd6tq+jyc45MoLeAfB0REtDqsqCngWFO4DocDRqMx4r22traIFb0WiwV6vT6h4zqdDi0tLeGf0P6ALS0tDH9LMF6xC1+c/DQK8jNjFTByChHIK0etfBDnBV+6qyEiIkrYigiADocDra2tMJlMsNlsaG1thdlsDh+3Wq0wmUwRbVpaWiAIAsxmM8xmM7q7uyPOWeh4SGjvQQAwGAzcEHoJhPz1+E3wJhTmrsqB5UWZ2HsAP5l6B84JE+kuhYiIKGEraiPo1YYbSUY6/dIf8a3fdeOzn/k7XFszfw/Gy5Eoitj81Udx/zuuwqf21KW7HCIiSgC/v1fICCBdHlSHH8QnFI9l1Aig7PRz+H7Oj3DOyRFAIiJaPRgASTo+NzzIz6gAiIkR3Dr1BEZHLqS7EiIiooQxAJJk5H433KIKhTmZshM0wlvByFzRnz5DRES0EjEAkmTkkx6MyVTIzcqgX6uSmfv+8kdPpbcOIiKiJGTQNzWl2kDBFpxS1mXWhsjZ+ZjIXYOywFl4J6fTXQ0REVFCGABJMg/VfgVP5DWmu4xlN3DtfXg2uA3n3XwmMBERrQ4MgCSNYBDBsSEU52TQ6N8lWbs/ieeDW3HOxQBIRESrAwMgSWNsEP/4+h3YIzuY7kqW3RrZCD6geBrnXdwKhoiIVgcGQJKGzw0AEHMzYwPo2bIGDuF7WT+Ea+hsukshIiJKCAMgSeNSAJTlZV4AROkVAICpoZNpLoSIiCgxDIAkDZ8AAFDkl6S5kDQoqUUQMiiFvnRXQkRElBAGQJLGpRHArPwMHAHMysVo9hoUjp9OdyVEREQJYQAkaWy7E3vEB5CnKkh3JWkxWPEWnPXlYjooprsUIiKiBTEAkiSCkOPcZB4KczPoMXCznN3zbfxo6t0YGvWnuxQiIqIFMQCSJPzPm/AvigdQmKtMdylpsa44DyXw4JxrPN2lEBERLYgBkCQRPPcytshPZ+wIYI3zGbyc+2mMDPA+QCIiWvkYAEkSQa8Aj6jK2BHA/MqNAADvheNproSIiGhhDIAkCdErwA0VijI0AEK9AdOQAyP2dFdCRES0IAZAkoTM54ZHzM/YKWAoszGirES2h3sBEhHRyscASJI4uvEe/H76hoydAgYAT14NCr18HBwREa18DIAkiaOljXhZthl5WYp0l5I2T1/7XXx28vPpLoOIiGhBDIC0dMEgqh0HcGX2CGQyWbqrSZuysgp4/EF4fIF0l0JERBQXAyAtnd+Ntx3/Jq7LPpPuStLqCvEUfpf9FQz2cyEIERGtbAyAtHReAQAQzC5KcyHpVV5SgmvldoyePZzuUoiIiOJiAKSl87kBAMHc4jQXkl6adVcgICowefFkukshIiKKiwGQls43MwIoy/AAKFdmYUBeCbmLU8BERLSyMQDS0inzcCRrK2T5mnRXknbD2dXIG+Xj4IiIaGVjAKSlW/8m/L3qO8hSlaS7krR7Zu3H8KvspnSXQUREFBcDIC1dwIcJrz+jN4EOCazbjb+O16a7DCIiorgYAGnpnjLiQf9nMvcxcLPU547hrokHMSlcSHcpREREMTEA0pKJXgFCMJ8jgACqCmT4gtICwWFLdylEREQxMQDSkk1NCPCIKo4AAiitqsOkqMD4hePpLoWIiCgmDtnQkk1POOGGiiOAAKo0hegXKxAc5l6ARES0cnEEkJZMnBDgFlUoYgBEbpYC5xRVyBIc6S6FiIgoJgZAWrJjup+gfepDnAK+5JDqLTiSfXW6yyAiIoqJQza0ZC5RBQGFnAK+5PCa9+El/xRuS3chREREMXAEkJZGFHHVE5/GHvmrHAG8pKZIifzhV8PPSCYiIlppGABpaaZ8WHvhLyiXuaHKVqS7mhWhPn8cJu8XIZ55Id2lEBERRcUASEvjFQAAk1lFkMlkaS5mZVCvqYVPzOJWMEREtGKtmJu2HA4HrFYrOjs70dXVlVAbo9EItVoNABAEAS0tLQkfFwQBBw4cAADY7XY4HA7s27cvfD4l6NI051R2cZoLWTnWlahwWqxE6eBxFKS7GCIioihWxAigzWaD1WqFIAhwOp0JtQmFO71eD71ej4aGBhgMhoSPt7a2YufOndDr9ejo6IBGo0Fzc7Pkn+2y55sZARRzitJcyMqxTp2H02IlxGF7ukshIiKKakWMADY0NKChoQEWiyXhNu3t7ejr6wu/1ul0aGxshMlkSuh4aMSxoaEBAFBfXx8eEaQkFFfj9xWfhR+V6a5kxVDnZ6FPVoOtgfMpv1a/cwIPvXwOE5PT8AVmfu7YXoW3XFGW8msTEdHqtSICYLIcDgcEQYg6XWu1WqHVauMe1+l086aZu7u7odPpUlbzZau4Gg+rPgCFmO5CVg6ZTIbO4k9gQFuOr6X4Wt959CishwexpjgXuUoFnBOTePmMgEf/7kbek0lERDGt2gAYjVqthiAICx6fy2KxQBAEdHZ2SlpnRrh4FJuEp3Gu8m3prmRFqSrJx3nBC4gikKIg5hyfxOOvD6D1nVfhnhu1AIC/Hr2IT/y0G6+f92DbOt6XSURE0a2IewClotFo4t5DGO242WyGw+EITwXH4/f74fF4In4y3pE/4B7hP7gJ9Bzri5T4tqMZOPRgyq7x0MvnAADvv25d+L0bN5ahojAHnT39KbsuERGtfpdVAFxoAUm043q9Hi0tLaivr0ddXV3UEcKQ9vZ2FBcXh39qamqWXPOq5xXggYqbQM9RVVYEr6iAePFISvoXRREHuvuh21yJ0oKc8PtKhRzvb1iH3x86D//UdEquTUREq9+qDIBarTbq+4IgQKvVLnhcEAS0trZGhD2dTgdBEGC1WmNet62tDW63O/zT389RFvgEuEUVRwDnqCtV4dh0NQIDqQmAh866cWxwFHt31QDjw8CUP3ysqaEawkQA/+/IxZRcm4iIVr9VGwDVanXUe/10Ot2Cxx0OB4xGY8SIYCgMxtsHMCcnB0VFRRE/mU70CXAF8zkCOEdtmQonxHUpGwHc392PtcW5uGljOdD1NeA764Gf3QE8acTG/HFsr1HD0ns2JdcmIqLVb0UFwFhTuKHANltbW1vEaJ3FYoFer0/oeENDA1paWiJGCvfv34+GhgauBE7SVEEVjgRrUMQRwAi1pSqcCFYjZ+ws4B+TtO+JySk8fOg87rpuDRRyGXDD54BbvgpkqYBn/wN4yIDmHdV44vgQLo76JL02ERFdHmSiKKZ9Aw+HwwGTyQSr1QqbzRa+Jy8U2MxmMzo6OmC3R26sO/tJH3a7HR0dHQkfFwQBZrM5/Dp0PJkngXg8HhQXF8PtdmfsaOCgx4c3ffsv+PHHd+LtV3EvwNkav/V73LVVhXve83ZALt3fWp09/bjf8gpe394JVYEauOPf3zh4aD9w6EG43/dL7DI+jS/dugn6m+oluzYR0eWA398rJACuVvwFAk4OCND9+7Po/PT12FWrSXc5K8qHzC9AU5CNH3x44RXmyWj+0XOoE/thHDQA7/4esPOTUc/73P/acHxwFI/93U3cE5CIaBZ+f6+wKWBafWp/vA2fVPyZi0CiqC1T4aYzPwRe2idZn6eGx9F9yoXPK38LFNcA1949/6RgEDhpxd7rKnB8cAyvnnNLdn0iIro8MADS4gWnoZwcxThyuQgkirqyfFR6T0I88bhkfT59chibFf1Yd+5R4KYvAcrs+Sc5HcAv78SewPNQZSvw7MkRya5PRESXBwZAWjzfzMgSt4GJrrZUhaPT6xAclG4l8IuOEbxDcxGo2AJc++HoJ5VdAWzYA3nvT3BNtRoH+12SXZ+IiC4PDIC0eL6ZrXM8UKEgmwFwrroyFU4E10Hh6ZdkJbAoinixzwnvlmbg088Aijijrjs/AZx+FreUu3CwP/bm5kRElJkYAGnxLo0ATmUVQi7nIoO5ajT5OIHqmRdDx5bcX9/wOErGTmJPlQKQK+KfvPkOIL8UuvE/Y9DjxwW3d8nXJyKiywcDIC1exVY8cN1vcDGnNt2VrEi5WQqMFV4B6/q/A4rWLrm/F/ucMGaZ8ebD/7Lwycoc4MYvQbNhKwDg4BmOAhIR0RsYAGnxlNnol61FTp4q3ZWsWFXlGliy3gMUVS25r8PHj+NauR1Zm29PrMH1n0XRjQasLc7lNDAREUVgAKTFO9GF2+3fQAEXgMRUW5aP3MFe4FXLkvoRRRF5p6wIQg5svDXxhhcO4SOaY3iZAZCIiGZhAKTFG3gFm0ef4wrgOGpLVdjieRqi9WtL6qff6cUu/4twl10HqEoTb9j7U3zEbcKrZ92Ymg4uqQYiIrp8MADS4nkFjMkKuAdgHNpyFY5OVUHmPgv4RxfdzwuOYQSgRO72O5NrWLsHJROnUBAYwfFBaZ9JTEREqxcDIC2eT4AHBRwBjKO2VIXj4tJXAr90yoX/LPsq8m78m+QabtgDALheeRSHznIamIiIZjAA0uJ5XRC4CXRcNZp89MmqIUIGXFz8htAD9oO4YUNB8g0LK4GyTXhH/gmuBCYiojAGQFq8bXfiIfFmFHEKOKYshRxlJWq8VvoOIK9kUX2cd43je96v4iMTv1xcEds/CIVmPVcCExFRGAMgLZq45X34zeRujgAuoLZUhf8ouh/Y/O5FtT/28lOokAkoa1hce9z4RXh23IfjF0cx5p9aXB9ERHRZYQCkRQsceRRrgwMMgAuoK1PhzLAHOH8QmBxPun3wyJ/hkRWgaOONi65hZ5GAagziFd4HSEREYACkJVA+9Ck0yntQmMMp4HhqS/MRdJ4GzG8FTj2bdPsNI0/BoX4LoFh80K7r+iQ+l/1HTgMTEREABkBarOkA5IFxeMBFIAupLVPh5HQ5plRrgVNPJ9V2cHgEE1PA9KbbllSDrHYP9mQd5UIQIiICwABIi+WdCRKCyH0AF1JXpgIgw0j5LuDUM0m1fa7fi/dMfgu1ez60tCJq92Dd1Fn0n+mDKIpL64uIiFY9BkBaHN9MAHSLKqjzGQDjWafOQ5ZCBrvqOuDCQcDnSbjtiddtuKqyAKWFuUsr4tJ+gFdMHMTFUf/S+iIiolWPAZAWJzgFT+EVGEERivMYAONRKuRYr8nHS+IWQKMF3GcTa+g+h5aTH8WnSl9ZehGFlfCtux4qmQ+HLyQeQImI6PLEAEiLU7EZv7vhNzgjr0Z+tiLd1ax422vU6BpQAff1ApVbEmrjsv0Wk6ICpVe/Q5Iasj/1ZzyivBVH0hwAz7omMM7taIiI0ooBkBZHFCFMBFCclw2ZTJbuala8N9VpcOSCBx5fABgbSqhN4NXf4TlxG3ZcVSdJDXK5DA0Vcpw4NyxJf4txdMCDW//tKew1PY9RXyBtdRARZToGQFqcl/bhU8/reP9fgnbXlSIoAv1P/gz4v5sWvg9wbAhlThteLXyrdFPsowP46VAzCs8+KU1/SXKOT+Ken/VgnToP/c4J3POzHvgC02mphYgo0zEA0uJ4XZiGHGre/5eQ2tJ8lBfm4BnvBkAMAmeej3u+6HRgAGVL3v4lQkElfDnlqBp7bdmDV2A6iM/+qhfeyWn85BO78MDHd+HQWQH3PfgypqaDy1oLERExANJieV0YlRVyAUiCZDIZdtdp8Pj5fKCwasH9AB15W3GD79+w/aorpCwC/rU7cJ3sOE4MjknXbwK+8fBh9Jxy4Yd370B1ST521Wrww4/swF+PXsSXf/sqt6YhIlpmDIC0OD4BHqhQzCnghL2pToNXzrkxtf4t8fcD9I/h0OuvQymXY1etRtIaVPXX4xqZA0fPjUjabzy/tZ3FL144jW+8dxt212kAnxuYnsLbrqrAvzZvh6X3LP5y5OKy1UNERAyAtFheF1xiAdR52emuZNXYXadBYFrEqcLrAOEMEPBFP/HIw3jfE+/EniqgIEfap6xkb3gzFLIgBk8dkbTfeH7+/Gm8/aoKfPhN64FzNuDfrwF++i7AK+B9163DpsoC/OnVC8tWDxERMQDSYr3re/iO+HEuAknCpoqZKfM/y98KfOkEkBVlc+cpP8QX/huvYBO2bayXvojqnfhS/SN4xi3tyGIsA24fDvYLePc1a4H+l4Cfvxco2QAMHQWe+i4A4Par16Lr8CD8U1wQQkS0XBgAaVHEoioc8ZUwACZBLpdhV20JXjgzBsgVwMUjwMCrkSd1fRXixaP4R/9HcX19aQqKUGBjVRmOnPcsy313XYcHoJTLcMtVlcCRh4HKbcDH/wjcYwXe/k8AgNu3VWLUP4VnTqRvexoiokzDAEiLMvXIF3GDeJCLQJK0u06D3tMuTAamgT9+Edh/98w9cQBw7FHgxR/hpY1/jxPyeuzYUJKSGhrHH8aDwftxwR1jClpCj70+iDdrS2fuFdV9Hfjob4GcQqBsI5CVBwy+jk0PvQu7ygL406sDKa+HiIhmMABS8kQRyoO/xAbZAANgknbXlcIXCOK1Cx7gvT8AJlzA7z8HiCKw/k3Ard/Ev7lvxo4NJcjNSs0TVirWVGGL/DTsfY6U9B/ingjgBccIvuP/JtDzE0Aunwl9s+WqAacDf1/yLLoOD2ByilvCEBEtBwZASl7AC9m0H25RBXU+F4EkY2tVEfKzFXipzwlo6oD3/QA48geg66tAXgl6192NF0+58LEbNqSshpJNbwEAeE48l7JrAMBfjg6iTuxH9dBTQJ46+knF64BtH8Au1x8x5pvEs3ZOAxMRLQcGQEqeTwAAuKHiRtBJylLIsWNDyUwABIDNdwC7DcCLPwImx/Hffz2JKyoKcOuWNSmrQaauwYi8DFkXelJ2DQB47PUB3Kd+DsgvA658V+wTGz6GrLFzaFKfwJ9e4WpgIqLlwABIyfO6AABusYCLQBZhd60G3aecmA5eWoRxWwfwhaM4PDyNvxy9iM/eXA+5PLXPVx4ouhprPK+krH/v5DSeP34Otwb+Clz7IUAZZ6S4eidQsQUfLn4Vjx8eRIBPBiEiSjkGQEpeXgler78XZ8VyFOYyACZrd50Go74p9Jy6NAookwGqUvz3EydRXZKHO7ZXpbyGvmvvxz0Tn4N3MjVbrzx5fAgN068hd8oNNHws/skyGXD3b5D1nu/B7Q3gOfvyfBAuwQAAIABJREFUbVJNRJSpGAApeUVVeHr9Z+DPLYMixSNVl6Pr1pdga1UR7vlZT3jrE8fQGP746gV8+q31yFKk/l/L9VdsxUVRjWODoynp//HXB3CubA9wn21mxe9CiqqwpUqNqzRyTgMTES0DBkBKnus0ii88zwUgi5StlGO/4Xo0bCjBx3/yEg709ONHT9pRXpCDph3Vy1LDpooCdGSZMdHzoOR9B6aDePLIWbxjSyVQmvhm1rKn/hW/FL+Mx1+/8Mb0OBERpQQDICXvyMN4//H7ef/fEhTkKPHAx3aieWcNWiyvwNJ7FvfeqE3Z1i9z5WYrcU32eRSc+YvkfXf3OXHP1K/xmb5L29skqmYXyrynoPW9jmMDqRmZJCKiGQyAlDyvC+PyAu4BuERKhRzffv82tN12FbatK555Vu4yulB0Lda6DyYX0hLwwolBNCufRv6G62bu70tU7U0Iqjfgw1lP4MU+3gdIRJRKDICUPJ+AURRwClgCMpkMhrfW4w+f2wNVjnJZrz1ZfT3Kg0OYdp6Stt9jj6EMAmQLLf6YSy6H/LqP4l2KF/CynfcBEhGl0vJ+48ThcDhgtVrR2dmJrq6uhNoYjUao1TMbzAqCgJaWlqSPA4DdbgcAmEymJX2GjOF1wY0CFOetmF8fWoSSLTcj+IoMw69aUXnzvZL06Z2cxhXOJyEUaqFesy35DjbfAfmT38PIqVchitdDlswIooT8U9MY801BLpNBLpchRylftul5IqLlsCK+wW02G3p6eiAIApxOZ0JtQuFOr9cDAKxWKwwGQzjELXS8tbUVHR0d4f4MBgMaGxsTDp8ZLacIfeIaqPM4Ariaba6rwecDf4PbFdtx2/9v787joq7zB46/Zob7HFEuFZXBE2+8yztQO8wyPGotu4S2e9tdyd22ayt+6FbWdoGVXbYpZFmWFpNpaZkC3mcypngrMCDXcMz39wdBotzMMDDzfj4e/DHzme+H95fvfJn3fE4L1bnjWA4DyKS8143Nq8C/D9tmp7HlvZ1kni+gZ4C3hSJrvLziMq575UdOGourn/N00fBJzBgGdvVt9XiEEMIa2kQXcEREBDExMeh0ukYfEx8fz+zZs6sfR0ZGkpSU1Khyo9FIRkYGRqOxujw2Nha9Xo/BYN39Ue3C9KX8s+wemQTSzvm4ObPPL4pfcj0tVufW33K5VfMf/Kb9o3kVqFQMDQ3EU13KVhutB/jG90fIKSzl9dsieGteBG/8KYIuHdx55st9KBYeLymEELbSJhLApjIYDBiNxuru3UtVJXH1lQOkpaXVSPaqks9Lk0JRu9JyM4WlFTIJxA6MCTIz5OCLYDxukfrSM08zIrQTajevZtfheX4XO1wWcOzQTovE1BRZOUUs3/IbsRN0XD8omGkDgrluYDBPTe9P2rFc1soahUIIO9FuE8DaaLVajEZjg+VarZbc3FwiIiKqy6oSw/paIU0mE/n5+TV+HJHm1UHcrvlWJoHYgT6dO3Jj4WoqMje2uK6S0nJeOL2AezVrW1ZRYDgqlRqvrA2t3uK2+JtDaD2ciRlf8//A1T07EdkvkP9bd5CSMuvsniKEEK2pXSaAdfHz86t3DGF95fHx8SQmJtbaanjpa3x9fat/QkJCWhxzu2M2o754kgo00gJoB/qGhrBP6c7FgxtbXNehXVvprjpLQK9hLavI2Z28oNGMKN3O8ZyiFsfVWDuO5/LlrlMsmhSMx/LJ8NbYPwo/nMlTYz04d7GEZT/IMBEhRPtnVwlgQxNI6iqPi4sjNja2esJIXRYtWkReXl71T1ZWVrNjbbdM+agUM3mKp4wBtAP9O/vwizkc56wtLV4PsGjPGvLxoNvQKS2Oy3PA9YxQHyLj0LEW19UYiqLw3FcH6BvkzYxTL0POURh1X1UhnDtAyMZHuWtMV97YmMnZ/JJWiUsIIaylXSaAdXXTGo1GdDpdg+WXSklJISwsrMHkD8DV1RUfH58aPw6npHKMpBFPtNIC2O55ujrxm/cwPEvOQO5vLaor6NR37Pccg9rZtcVxuYdPQ61SOHNoa4vraox1e8+QfiyXV8J/Rb03GW54GYbOqyxUqWDWe3BiO4+5f427i4bF6w+1SlyNpSgKecVltg7DriiKguF8AdkFJluHUq+yCjPlFWZbhyHaoXabAGq12lrH+kVGRjZYXqVq3F9V8lff+EHxu+LKBDBP8cRHEkC7YO42mk88bgOn5iduJXnnCCrLojB0qmWC6tCdJYPXseJcD8vU14B3Nh9lQpg3ffYsgUFzYGB0zRd0GwVjH8NtyxKeGmZizc6T5BaWtkpsDSkwlfPoyp0MefZblnxzkLI2nAxcKDBRaCpvs7OpK8wKPxw+z9Nf7GP8ku+Z/OImJizZyPItR9vc/tR5RWUs1R9m2L9TGfpsKgs+SOO9LUc5cq7A1qGJdqJNrANYpa4uWoPBQEpKSo2FnBctWoRer69O3lJSUmq04jVUnpGRQUZGBtHR0dWzhleuXMmiRYuscWr2o1Nv1oxcwcmfK2RhXDvRu1sXntx3IzM9gmjutJ5d2RrmmxL5dOTYhl/cSEN69SDxl+2czC2iSwcPi9V7uTN5JaQfy+XFWYOhx1fg2an2F058HAwbmeydhVnpxvp9Z7h1ZOtu33e5A6fzeWBFBmfzS5g7IoS3NhnYfCSbV+cOoXtHyy3v01JVXezvbD4KgLNGha+7C1HhgTx/0wDUatss+H0pRVF4bNVO1uw8RWdfNyb3C2CSzosfjuTy7Nr9fLbjJC/cPJABXWy7FmRecRnLfjDw/k+/UVZRTly4ER9XDV+ed+L5r89RVrGf2PE64qb1bRN/16LScj78+Rin80pQq1Ro1ODu4sScESF00brbOjyHplLawFcxg8FAYmIier2ejIwMFi5cWKNbNikpiYSEhOodO6pcutNHZmZmjYWd6ys3Go2EhobWuuRLU/4c+fn5+Pr6kpeX51DdwS99e4hVaSfY+o9rbB2KsID0Y7k88tYaVk3Ko3PUw03bv/d3r6UeIHHLcXY+OQWNhT50jFn7KX97Knsmvs2kSRZqWazFe1uOsmndSpb+LRZfbYf6X1xRDhon5r39C2ZF4eMFo60WV0NWbj/Ok2v2ofP34vXbhqLz92LH8Vwe+WQn2QUmEqIHccOgzjaL71JvbswkYf1BHo3sRWgnT/KKyzhpLGbZDwbuujqUf90QbusQ+e93v7I09QDLIxXGXXNj5S40/x0OJXmcCF/AA4eHsOdcGS/NHsJNQ7vYJMai0nLmJm3l7NkzvNBjJxMKvsIp9/deqy7DKZr/DR9tPcab67YzYXAfFkcPxsXJNh19iqLwxa5TxH99EM+i4wzyNVGuqKhQ4ECxllPlPsSO1xE7IazVt8EEx/38vlSbaAHU6XRXJG+XiomJqXWM3uVbuzW2vGoZGNEMhk2M+fUjvvW43daRCAsJD/ahr/oEnX9aAsOng1/jF2QH4OIZ7v5pMqbApyyW/AFoO/eiUFVO+cFvwIoJ4Ld7jvOWyxv4pJfCNU/W/2KNE5gK+LvXOv68S8e5/CEE+LhZLba6HDidT9yne5gzPIRnZvTH7Uw6fLmCocPu4utHxrEwZRd/S97FiB5+BNogvkutSssiYf1BFo31JdZ9PZQ4w+A54NGDLlp3nlyzj2BfN+4d18T3nQV9tfs0r6XuRR+8jNDNW2DkAfDpXDkWdPdKuqYn8LmbL2u7RvPP1eX0CfKmX3DrJg3lFWYe+ngHR84V8PlcHb1X3w/hM+Cm18FNC6UFeLg4EdPlGHd7/oUn9t7BPQUm3rx9OF6tnGAdPnuRRav3kH4sl3cCU7imdDUU/lFuuuZ5XikczIof9pCy/SiP3zCIGwe3jS8rjqRNJICiHTmZzqDs9fgG3m3rSISFuLtoyO0YgTlfjfroj01OAMt3foLGXEbH3hZuDdM4c9R3JF3O/2jZei9x7mIJAVnr8XHOg8G3NvIohUGGZcxziuTrPeO48+pQq8VXl9e/P0LXDu48N6Mfzs4acPOFHR9B+vt4RdxOwrTHGZeZzSvf/coLNw9s9fiqfHfgLItW7+HhQQoxv8ZC4QVQqaH3VHDTckf2K7gMHszjX0GQr5tNWix3nzDyj+RtfNbxDXpc3A23rgTv4MrC0HGVP+P/jmrzy1x/ajsf+Y3nzx+l88VDY/Fxa51x0Iqi8MwXewg58hHzb32M3v27Q9ivldf9cl1H4jRwJv+34y3WZe3m3rce4e37IlstCdx3Ko/4pI/w9tTy8b1TuapCDcVR0GUYKGZQzLh6BbLQsxP3F75O6YH1vLTqeo5f+DMPTO5ts/2/HVG7nAQibKg4lwK1lywBY2fCunXhgFNf2L+maQcqCqa0FaSahzEq3PItOGVhU+hb8Svnz56weN0A3+w7y+0aPWXdxkOnXo07yNUb1ZA/Mc95I+t2tc4yNZfKPF/AV3tO89+QTTi/PQnKTeDfB/55Fq5dDPu/wDtxFC/3OcDK7VkYzttmUsBJYzEPfJzBPaG5/CXrIVQuXvBQBvzzdOWXjIpSyDvJnF//yr9D9/HYyl2kH2vdnpnsAhMPvLeF991epp9pD6rbVkGfaVcOg+jQHaYvRR2zgcV3TILCc8StTG+1ySyJ3x9mUPoTPOX0PuOdDlQ+WVvyB+DqBTNeh+jlTHHdy39yHuS5D79ulUksB8/ks3JZAst5krcGHOSqnp2gz7Uw5LbK92hAPwjsXz3O1mvsfXQIn8hzzsvpuzGWl7/8pc1OELJHkgCKpikxkocXWnfZBcSeDOyqZUXJ1SiG7yH/VOMPPL0Tz7zDfOcaSd8gb4vH1XXEdNQqheMZqRavG2B/+maGqQ/jPHpB0w4ccS8+ZiPBJ77hpLHYKrHV5Y3vM5npuZuhh5dC6Hgw/74zicYJRsXAwztg7COMiYomwNuVF1MPt2p8f8R5BHdnDY9ePxRVyCi4ax34XjJ2ztkN5q5ANWgOt59+nse1ev71+d5WnW27VP8rFeUmwgPdUc1LAd2E+g9QqeiudWWtTwLXH3mSZRut/7fduO8EPb5/gFuctqCauawyQW2MATPR/Hkzzt1H8XWmif98a92li349beSXxAd5Vnkd86C5uEX9q+GDAvtXntNtyYx1MxCdNo8nU7a3uRnX9koSQNE0xbkYFQ98pQXQrgzq4ssX5aPI6T0XzOWNPzBrGxdUfrj3vcYqXTf+wd2Y4/E2a0pHWLzu7AITe0/m8lvwtdDnuiYG1pvyHhO50/lbvtrdhIS5hbJyili38yhPO38IYZNh6vPgctkMaQ8/GP933Py6EDc+EP3uY+w5kddqMQKcMhbzfdpuYsd2xaNzONz2SWVcl9M4V7ZWjX2MuwvfJuzselZnWKe193JHzl0keZuBOycPxuWer6FHI2ewa5zwuvZZrtOk0WnDY+w6Xv8GBC1hLDSRn/IQ12h2orp1BQya1bQKOnQn8O4VPHhtBF9t+gn9Dz9YJc6j5ws4ljSXecqXFE1+DtebXwOnJjQS9J6C2wNbOBfxCCsyzvOXT3ZIEtgKJAEUTdN7Gl+bx8g2cHamb7A3Jo0nX/eIA23jlzbJ6nU744v/w/i+wVaLrYeuD78Yclq8U8nlvt1/ln1KKF5/er+y9ayJnK55gk1dYvhyZ+slgG9uyuQR17V4mc7BtUvqn7FdbmLGtrk84/MFCesPtlqMAG9tPMIrzq9z78mnGn6xSgWRT8Eda6D/TSz55hBFpU34EtJMq1an8L3rY8zvqzR95nu/G1BmJjFTs5ndH/+T0nLrrL341Bf7OKYEUDj1JVR9rm12PQvG6Xir4yqGfDePA7ssu7h6bmEpd723na2uYyi45WM8xj/UrJUE0IYw/KaHeO22CLrtf4utr86noszCa22Wm6CsGCpk0XSQBFA0kXnwn/igZJyMAbQzrk4aBnfVsuXwaUhbDmf2NHxQYTabDp6mVO3O1T3rWDvPAiYEl/F6biwXD1u29eLUts9YEHyUTl7NXAA7ZAS9xtzInlP5HL1Q2PDrW+h0XjEpaSfoPHAiqqnPQ6ee9R/g5Ipq+N3MKf2cgsytbP71gtVjhMp1FS+mrWQE+3AaE9v4A3UTWXhtf7oX7+ft7w9YLT6Anw4c57ZT8bh16IJrp+ZN4tEMiubc8L8yp3glH63fbOEIQZ9+gDW7TtN1xr/QjrmjRXWpVCrCYj4g37kjwZ9Fc/5ImkViNBlP8/lb/+RicRl3LPg7vgObn6RWuW5gMFEj+jMydy2GpVMxF7a8hdVcnMfe1+bCcwHwfBBvP3U7A5/+psX1tneSAIomKcrcjL+SI2MA7VBUeCCbjuRg/mEJbH+74QO+eJABP97P8B4d8LbibMjB4X3Rqgo4v2OtxerMLTAx4/xb3OqyqUX1TO6Ux4euCaSm7bdQZHVb9sNR3J3VTLxuDoxqZGJ11cPQeQj/9Uhi+Q+t0wq4fMNuFjl9RFmf6dAzsuEDLhHiWszHLs/hsyWeM3nW2W+5wqxw4bPHCVYb0d72Dqibv6B9wHVP8PHgD4n/uYDDZy9aLMbcY3sZ8WUkj3c/zE1DLLPmoKtPAD4x6zitCsB1xU2UHM9oUX1K3gmMb0Rxff4q3p0dSreOllusffCMR9g69l06Fhwm55WxmE/tanZdpWf2k/PSaLqf38iGbg/z0+B4uk6Yz5NtYO1JW5MEUDSJx6rZ3KD5WVoA7VBUeCDFZXC0y42wd3VlV0ldCs6j/PotXxb2ZWKfAKvG1aWDJ9s1EXge+85ide7avJaeqlP4jv9zi+px8/FntHo/zrs+slBktTOVV3A+fQ1rfJbgrWpCYqRxQnXTG3RWzjLo6Dv8ZuWWyrP5JQTuWEoHdTHO18Y3vQLPjlRMfII71V+xZvUKywcI/PhNMjeWfsW50f9A1VArakPUauZOn0aPDq5s/eAJKgpbPotZKbxA6UezuYCWWbNvt+jY2k4BQXDH5xw2d+Gdb7Y1f8btmT3kvxFFmamYvVM/YXDfRs6gb4JxUTexNfJTzpU4kfHxUxSamjgsQFHILynjvs9Ps83UjfRpa5h897+56ub7mRZ1LbOGh1g85vZGEkDReKVFqMuKyFF8ZAygHdL5exHm70ly2Tgw5cPBr+p+8Z5kFFR8WjqGiX38rR7bueAJBBZnQp5lJgi4713BCU1XOvSb1LKKPDtypus0JhSu52RukUViq83Gg+d40PwRAd6u4OLVtIMD+mGetoSdzkNY8Yt1l61J3GRgpzqc8inPg7Z5H7BuYx/kdMdRzPjt3xw0WDbeotJyPko7x8++1xEy5WGL1OnmrOHFawO5sSCZ0+/MhfIWjFsrLSJ72c2oSy9yYtpyOna0/NCKfqHdOH/LZyzJDOH17w7CsZ+bVsGJNMqTruFEsTPfjX6fyVeNsXiMVa4bN5qs6LU8mH87t7z5E+cy1sKJBrqvy0shbTmmpRHc/8aXpJ+toNNd/2PimFFWi7O9kgRQNF5R5RiibCQBtFdR4UGsOuqCEjIGdtbRApN3ErYs5YDvONx9/ekTaPnlXy7n3X8K5Yqa4gMtH7dTcjGHwRd/4Hj3mc0brH4Zv6vvQqc+w86frDem6OBPX9BHfQKPyMebFbPzyLvoNWIKq9JOUFJWYYUIIb+kjE+2H6fHVbNwH31v8ytSq/Gf9w5eqlKOr25gZ5YmWr75KJtKdHSd/w6oLffxN6j/ANb0+T/8s9PJWzH/j6V5mih37ZN45B5iZa+XmDB6pMXiu9y1gzrzWFRvjmz4AJZPgzUPQvGVW6PWoCgoisJLe9x40XQzG67+iDumXW21GKtMHdSNDx6YQklZBYe++A+8fQ28MQb0z0DWdjD/PgHn+C+weSnm/w5DWfsXUnMCMZvNpNw3hpGhtcxAF5IAiiYo/D0BVHykC9hORYUHklNYyuE+sRBRx8DzEiN0COVJ0x1M7OPfKiv3R/QO5brSeLb6NHG5llpsO2ZkSfkcAsfOt0Bk4Nl7Euc1gbjtX2WR+i6XV1zGkBMfc8GrT+OXKqnF/D7lvFH+NKm/NH88VX22fLuaV1jC7UNb/mHr1CGEPRMS+duF6/gp0zKTVy4e/J4Rm+Zz77AOhPhZbrxalTmz/sRinzg8j66n9POHmjxrvaSsgnt/m8w/PJ/hntkzLR7f5R6a3BPP4bfyeNm9mHZ9ivLGaDjw5ZUtmHknYfPLKG9exTsrU3h10zG0Uxby0NRBrbZrR+9Ab9Y8MJZ3Q+K5r/RRNuQFU/LLu/BOJGQfobTcTN5P71Kx4QU2FnTlhvLFHJv8Gu89OpNerfAFtb2SBFA0XrmJQtcAjCrfVt9bUrSOoSFaOnm5kpLXB/rfXLlswr7PKwvN5srHgf05PmM16TkuTOht3fF/Vbp39MDo1ZNfjrV8Pbtvfi0k1TcaXWiYBSID1Gq2jXyFv+beQl6x5ZeX2JBxgAjVIZyvfqBFLZZdOocw2Pk4qs0vWjC6SuYKMyE7X6K3RwH+nSzTbTl64g2Edu3Csq+2oBS1cCZoWTFlnz2IE+XcEzXEIvFdzs1Zw/w7H+Bf3E+G4SxKY9fTNFfAlld587Pv2JOj4r47/oS7S/MnpjSWSqXiuZsHMmD6w0SWLGZnaVdYOQ9O7ah8wa5P4P0b4eX+mL+PJ60oiM92nePfM/oTO8FC904T+Ho48/Zdo5k+9z5WBC9iYOHrzC1/ipnv7qXfk+sZuXMaYUXvsrzz07z+6J94YFJPXJwkxamPfIqLxus+huWjvqZky2+yX6OdUqtVRIUHkLr/LP+4rh+qA1/Cp/fAkdsrt546tQPmf8nGX8/jpFZxdc+OrRKXSqViYncXpmfcA73/XbkIcjMo5w4wYM//4T3wIYu+hyNGT8T4/QY2HjzLjKFdLVYvwMp9hazrsoKkEc1v/QPAXcuZ/guYsvsV9u/fS3j4AMsECOz9cTWDzAc5PPY9i3SrQ+U1j5vak5AP7+TU/0bT5Z7mT7S5uP5ZPEvOsmfYx0R4u1skvtp06+jBxFkPMffDNJ78+QR3B2VC0CDwquOLUu4x+CwW5fgvZJfN5183PErfIB+rxXc5lUrFvNHd6R14PQs+CmSQ0xG6bVfT58xxJp85hGtxOav9/sJLp8Lx1HTgH3P6McNCs5KbQ6NWcf2gYK4fFMyFAhNf7OzPidxibgnwpKe/Fz0DvOjY3GWdHJAkgKJJjEVlsguInYsKD+R/27I4cq6AXgOjK/dsXfsXKC+BKc+jqNSs3X3a6su/XG5gWDc8f71A+Y7/4dTMBDD7x3eYUrGJw+FLLBpbsK87T/p9Rw/96zD0a4vVe+rcefYZTvDUrDHg1PIPttDrHyN/z9vkfxsP4RaaZasoeP+UwH5NP/pdNcMydf7uql5BJHW6k5isJVTsX4sm/IamV3IyHc/0t3hdcyt3TW3e+6YppvYP4t6xOhK+3sttHR7HrSwPJvwdRt1X8xruToavHqNQ7cU95f+iQ78JzBvd3erx1WZkqB9fPDSOpfpAfj6exwdpezArw4Hh9A3y5plZOqYP7tymWtQ6ebly99jmreEoKkkCKBpvw3Pcuv9b0r2asbyDaDeuCuuEh4uGb/efrRw/M+Q2CBoIJzMg4g6+2XeWbUdzePuO4a0a12hdR1aVT+CvB9ZA8RJw1zatgooyPA6mkKIaz609gyweX0jXEAYffgfThd9w7dTDInUeX/cKG13fx6WPZdYZ1Lh5sy/sXvoceYfc3Fw6dOjQ4jpPHdlJsOkoP49+k3ALTqyoctXMB9En6hn7+YNoAvtBx6Z1P546nskpc098Ix9rtaErcdf2Zf/pfK46+hTv9tAzWP8Mqm3LoPMQmPNR5YSKtX/hkO/VzMqKZtqw3rwwc6BNe1Y6a91ZHD0YgOLSCg6cyae8QmFEjw7S42On2k46L9q+vJMoFWVoZQawXXNz1jChtz/f7j/7x5NBA2HYfIrLzPx77X4m9fHnmn6tM/6vSs8AL753j0RVUQp7U5peweFv8CjL5USPW3DWWP5fX7dxcylUXDm16R3LVFhRRthvH3NYOxZvL8sNZO87/VGmVbzE/3ZlW6S+5YddmaZ6k9HXWGfiwoCuWr7r+zRnSt0pXzm/8ZMrKsqoqDDz2J4Q/uYZz9zROqvEVxtnjZoP7h7JvElDmHl0Bo/5vU5ByETo1Jui0nKO5RTx3y6LmZZ1BzFRQ1gcPcgq78nmcnfRENGtAyND/ST5s2Nt5x0n2r7C8+Tii9ZDdgGxd1HhgezKMnLgdH6N59/clMn5iyaenN6/1T8YVCoVwwb04yf1MJS9q5t8fMn299ll1tF/yGgrRAe9uwax0WksvoeS/1iaogVObvkYfyUb9ZiWLVZ9uU5aXyKH9OaLLbsoPdOy3UFKTu3ni+2/MmXkAKtOXPjL9NE8oH6C55wepFHpn7kCVs1n93sP88vRHF6IHtLq3ZdOGjWPRfVmVewYthUGMmzn9fTfNJLwJ79hwpKNvHa4A0vnDOGha3pJkiVsQrqAReMVXeC82V/WAHQAU/oH0SfQQPSbP/HSnCFM7R/E8ewi3tqUyYLxoYR28rRJXNcOCOZvW+9g2aQpNHUKgz44hg8O7WGZlWYuq1QqLvS8Be2hDZjP7kMdPLD5lSkKmq2vsZVBDGvp5I9a3DMulOxdsRSsdMLv4R+aN3GjtIjyD29hUUUoEaOa0SLbBAE+bvwlOpJ7P0ij55bDzFPWwpiHwKmOL6Pf/APl8DpeL/0rsePDuCrMentVN2R4Dz/WPTqOFVuP46xR4e/tir+XKzp/L4J83WwWlxDSAigar/ACZyu8ZA1AB+Dl6sTq+69iXC9/Yj9M5+XUwzy7dj8dPV14YFILt89qgVGhfpR4BPHVQSNUNGHJFUUh5YQvmu5XWXUSU9+RU5lkepFdZS2bCVxWZORYsSuHe91jla7B3oHebAy+C7/c3Sj7PmtWHeYfX8Kl+BzpPWIsug9sXSLDA7ljTHc+X/8bs7L4AAAbpUlEQVQNyvfx8O4U2Pm/mlsWHvkO3rsBfnmLpS4xnA2ayGNRva0eW0N83Jz588Qw7h2nY8aQLlzVs5Mkf8LmJAEUjWaes4L3TBPx85QuYEfg6erEm/Mi+NuU3ry64Vf0B87yz+v74eFiu44DJ42aKeFBmHYko7w8oP79iqtcOELFf4dzNnO31cctDuvhR5FXd77YkQUl+Q0fUIdNx0qZU7yI4ZOstyDwuKiZpFZEYFr/ZOX6jk2RnYmyZSmJ5TcwZ1oLt9Nrgn9c14+LHYcQ5/40Fa6+8Pl98FI/yPiw8gXHt0J5CR+FPEtS0SRemdv6Xb9CtBdyZ4hGy/Xpy7GKTpX7kQqHoFKpeHByL969cwQPX9OL6wcG2zokpg0MYkN+MKqCM3BgbcMHbHyBkqKLHK3oyPWDrBu/k0bNzKFdmLHjXiq+/VfzKjl/iP0/riY8yJvwztZbE+7qnh1Zpb0X54KTsG1Zk441r4vjvKLlQNgCBnTxtVKEV3Jz1vDqrUP5PC+M20sXcWjWDzD0dvDpDMDpoQ+zyO8lnvi1J0/fGI7Ov4n7JgvhQCQBFI1z8SyqdQvpqjqHv7d0XTiaSX0CeCyqd5sYrH51WCeyXUM44TMUMt6v/8Vn9sLeT1mmimZ8vxCCfa23CHCVW4Z1JbVscGX3ZGHTtzEzbUhg5qn/MHuYdZNVlUrF1IkTeLrsDo5rm7akz8/+c/ib6V7+HNWCcY7N1CfIm8Tbh3Emv4SpH55gwZkZbFUP4dkv9zPhxc2s33uGJ28IZ/bwkFaPTYj2RBJA0Th5Wfjtew8vSqQFUNiUi5OaqH6BrCibCL/9CJnf1/3i75+nxLs7r+WO4vYxrbPIbu9Ab3YG3ky5Amx/u2kHG4/jfPBzlpuv58YI68c7fXAw6z2m8+ZBj8ou64YS1nITFRVm/rXXH5fe1zCwa+u1/l1qUp8Avn10PC/OGsyhMxeZm7SV5LQsHpjYkx8WTuLusaFt4suKEG2ZJICicQrPA5CteBPgIwmgsK1pA4JIyo2gsNskOPZT7S8qyYNzB1jl+Se6dvTh6lacCTp1RDirysdj/mVZ48YpVvn5DQrw4HzP2a0y1tbVScOdV/Xg0/STFKy4HT6eDaVFtb+44By8O41zy27BcL6AR67pZfX46uOkUXPLsK5899cJfHjPSH5YOIlHInu16u40QrRnkgCKxvm9ZUBx98PVyfoblQtRn/G9/XFzceH9bvEw+Z+VT146K7jgPLj5cuHOzTyf1Z95o7ujVrdei9D0QZ1533w9+Sqvyv1eG6PwAhXp77O8LJIbR7RecnXP2FB0/p7E5c5AOXcQVi+oXEfvUhd+hbcjUfJP8szF6UzsE8DgkCbuxGIlzho143r500EmpwnRJJIAisYpPE+Rxhs/H9us/ybEpdycNUzqG8BX+3/vsty7GhInQNY2WDkPXhsGRTmszDgDKg3Rw1q2LEtTdfB0oXf4IOY6v4ri36dxB7lp2RQwjy9cb2RCH3/rBnjpr3XWsHTuEFKNwazs/jQc+hqWX/dH13rWdngnCpzdWdY7ifXZgTxs49Y/IUTLSQIoGid4EBt8biZAJoCINuK6gcHsO5XPL4ZsCBoEF09XJiqndsG0BCpcfPj4l+NMH9zZJrvXRA/rysGzBfx6cDf88J/6tzArzKZUUfO3s1OYHNG31bcF6xvkQ9y0vjy+twsHrl5auc+y8nsr4Ml0CBzAB/0SeeGnQhZd25eIbi3fQ1gIYVuyE4honJ6RLHfxoLtMABFtxJTwQEaF+vHwJzv46uFxdJr/ZWWyMnguOLmy8cBZThqLuX1060z+uNz4Xv508nJl39Zv6X3sOXD1gVExV77w1E5Yfi0bBr5MbpE3c0bYZvbqXVf14PuD57hzmyvrH/nwjy7V/jezShXFk58d5P6JYcROCLNJfEIIy5IWQNE4p3agyfuNAB9pARRtg5NGzau3DqXCrPDoJzupCOgPw+aDkys5haW8uuEIA7v42mysmpNGzcyILjybNZjSkffD+rjKnSouZSqAlLsp6xDGP9K9uHVkN3oGeNskXrVaxYuzB2MqN3PtKz9y/4p0Xv/+CG/vLOTxzw8yb3Q3/j61kd3ZQog2TxJA0SjKFw9xc9GnsgSMaFMCfdx4Ze5QtmRe4LUNRwDYdPg8U5f+wPHsQp64vp9N47t9dHdM5Wb+VTgLekZC8l2QnVlZuH9N5XjFi2d4yedxzBoX/j7FtglWoI8bH90zihsGBZNdUMqbGzN57qsDTB/cmWdvHCBLqwhhR6QLWDSKUnCBc+ZehMkSMKKNubpnJx69pjdLvzvM4XMX+Wr3acb16sR/Zg0m0MYt1iF+Hjw9vT8LP91N5JwXiDLO+6PwyHeQdwLD2CW8uU7Fczf1aRMzWQd08a3e3cNsVjhfYCLA21WSPyHsjCSAomGKgqroAtmKN2NkEohogx6c3JO0Yzmk7j/LU9PDmT+mR6su+1KfWcO78t3Bs/z9y6N8e/9XBHj9/iXqxlepMCs88voW+ndWuHVkN9sGWgu1WmXzJFoIYR3SBSwaZspHZS4jW/ElUFoARRukUatYdsdwNi+cxF1Xh7aZ5A8qt1yLnzkIF42av67JxOxcuT+toiis+OUYe07m8eyM/mjaUMxCCPsnLYCiYaaLFLkHcc6klWVgRJvl5qzBzbltLlLu5+nCklmDmf/uNhZ8kEZ+SRkHT1/koqmcWyK6Mqy7n61DFEI4GEkARcN8u/LB6K85uOEI7i5t8wNWiLZuQm9/Hprck2/2naFPkA8T+wTQL9ibcb1ab9FnIYSoIgmgaJRz+Sb8pftXiBb565Q+/NXGM32FEAJkDKBojB0reGDXTQR6ySbrQgghhD2QBFA07OIpnMwlBPh62DoSIYQQQliAJICiYYUXyMVHFoEWQggh7ESbGQNoMBjQ6/UkJyeTmpraqGMWL16MVlu5zZPRaGThwoVNKm/O73RIhec5W+EtM4CFEEIIO9EmWgAzMjLQ6/UYjUZycnIadUxVchcTE0NMTAwRERHExsY2urw5v9NRlRec57zZmwCZBCKEEELYBZWiKIqtg6iSkpJCfHw86enpDb62Q4cOHD16tLqFDyoXXK06nYbKm/M7L5efn4+vry95eXn4+Pg0+fj2IuvwTu5cvo3n7r2FMWEdbR2OEEII0SKO8vldnzbRAthUBoMBo9FYI7mrotfrGywXTXPSKYRMpYu0AAohhBB2os2MAWwKg8FQ6/NarRaj0dhgeXOZTCZMJlP14/z8/GbX1W6YzQRufoKBqnACvKfYOhohhBBCWEC7bAGsi5+fX73j+Roqb0h8fDy+vr7VPyEhIc2uq90oMRJq+JhQpxy8XNvl9wUhhBBCXMauEsCGkruWTvZYtGgReXl51T9ZWVktqq9dKLwAgOLZEZVKNqsXQggh7EG7bNLR6XS1Pm80GtHpdA2WN5erqyuurg42Dq7wPAAarwAbByKEEEIIS2mXLYA6nQ6tVlvrWL/IyMgGy0UTFFW2ADr7SgIohBBC2Is2lQDW1UVrMBhYvHhxjecWLVpUY0ZvSkoKMTExjS5v6HeK3/mG8Lnz9Xj5yvIvQgghhL1oE+sAGgwGEhMT0ev1ZGRksHDhQsLCwqoTtqSkJBISEsjMzKxx3KU7fWRmZpKQkNDo8oZ+Z2M4yjpCA5/+hvsn9uTPE8NsHYoQQgjRYo7y+V2fNpEAtleO8AYynd7Pda9u5oFZ1zEzoqutwxFCCCFazBE+vxvSLieBiNZToX+Op52OofKeaetQhBBCCGEhbWoMoGh7FONxTij+sguIEEIIYUckART1cr6YVZkAeksCKIQQQtgLSQBF3UwXcSk1clodgK+7s62jEUIIIYSFyBhAUbfCC+S5BFOo6iq7gAghhBB2RBJAUTe/UJ7tuZLzFwpsHYkQQgghLEi6gEXdFIXfsgvp5udh60iEEEIIYUGSAIo6Kd/8k8fP/o0wfy9bhyKEEEIIC5IEUNSp9EImhRVO9AyQBFAIIYSwJ5IAijqVZ/9GluIvCaAQQghhZyQBFLVTFFzysziFP907eto6GiGEEEJYkCSAonbFuThXFFLiFYKLk7xNhBBCCHsiy8CI2rn68HhQEkWugbaORAghhBAWJk07onYaJ37I7USX4CBbRyKEEEIIC5MEUNSqZNdq7ip8m56yBIwQQghhd6QLWNSq8OD3jFXvwSQzgIUQQgi7Iy2AolZl2Uc5ofgT5i8zgIUQQgh7IwmgqJXzxRPkOAfj7eZs61CEEEIIYWGSAIorKQreJaco9e5q60iEEEIIYQUyBlBcqaKM913mUh402taRCCGEEMIKpAVQXKFM5UTCxWl49YiwdShCCCGEsAJJAMUVThv2MIE0mQAihBBC2ClJAMUVSnZ9zovOb9EzwNvWoQghhBDCCiQBFFcoyz7KKVUA/t6utg5FCCGEEFYgCaC4glN+FkbXzqhUKluHIoQQQggrkARQXMG7+BSlXrIEjBBCCGGvJAEUNShmM8cqOlIR0N/WoQghhBDCSiQBFDWcvVjKraZFKAPn2joUIYQQQliJJICihszT2ThRTs8AL1uHIoQQQggrkQRQ1FCS/jEHXO+iq+R/QgghhN2SBFDU4HpsIyfdeuHkJotACyGEEPZKEkBR7VR2PoNMGZSEXmPrUIQQQghhRZIAimq7f/4WH1URXUfdZOtQhBBCCGFFkgCKamePpJOj6YhX9+G2DkUIIYQQViQJoAAgr6iMZ8+NY/3kdaCWt4UQQghhz+STXgCwcX8WZnMF1wzsbutQhBBCCGFlkgAKAIq3Lme7x6MEempsHYoQQgghrEwSQEFJWQWdz/9AiU8P0DjbOhwhhBBCWJkkgIKfD2Yxkn24hV9r61CEEEII0QqcbB1AFYPBgF6vJzk5mdTU1EYds3jxYrRaLQBGo5GFCxdatNxR/LZ9HZNUZbgNvdHWoQghhBCiFbSJFsCMjAz0ej1Go5GcnJxGHVOVvMXExBATE0NERASxsbEWK3cUJWUVmLO2kevaBTr1snU4QgghhGgFKkVRFFsHUSUlJYX4+HjS09MbfG2HDh04evRodQsegEqloup0WlreGPn5+fj6+pKXl4ePj0+jj2srsgtMLPggjf2n8/nynv706tHN1iEJIYQQVtfeP78toU20ADaVwWDAaDTWSN6q6PX6Fpc7gqPnC1jy31eZd/5lVi4YLcmfEEII4UDazBjApjAYDLU+r9VqMRqNLS6vi8lkwmQyVT/Oz88HYM+mFLw8Paqfzw0YhdnJA6/c/bgWn61RR5FXD4p9QnEuycEnZ1eNsgqNO8bA0QB0OLMFtbm0RnlepwjKXXzxyDuCe2FWjbISjy4UanujKb2I9kJajTJF5URO8LjKczy3jbKLF2D72/wfeynpejVuAbL0ixBCCOFI2mUCWBc/Pz9ycnJqbdlrSnld4uPjeeaZZ654fuBPj+Ljqqp+PMH0EseUIF52fp2bNVtqvHZp+UyWlkczXr2LD1wSapT9Zg5kZunLAKS7PkxH1cUa5TNNT5Oh9OYJpw+512ldjbIPyqN4svwu+quO8pXrP2uUXVTcucn0DgDfusTRW32SE07dKJixAq8B14NKhRBCCCEch10lgA1NIGlp+aJFi3jssceqH+fn5xMSEkL23T9T6u1d/Xyyux+onVCVDOdChalGHX9y9uQ2F08oH8sF07waZZ4qNds8OgFgLvyJC5f9/rfctKBxQWUayYXy4hpl1zq5M83VGyrKuFBSczavAmzzDABAVZTKBaWCLgFdUKml5U8IIYRwRO0yAdTpdLU+bzQa0el0LS6vi6urK66urlc83zEwpPZBpN5uddYFboBv3cXe9YzJa7Be77qLvbvUc6wQQgghHEG7nASi0+nQarW1juWLjIxscbkQQgghhD1rUwlgXV2wBoOBxYsX13hu0aJFNWbspqSkEBMTY7FyIYQQQgh71SbWATQYDCQmJqLX68nIyGDhwoWEhYVVJ2RJSUkkJCSQmZlZ47hLd/LIzMwkISHBouUNkXWEhBBCiPZHPr/bSALYXskbSAghhGh/5PO7jXUBCyGEEEII65MEUAghhBDCwUgCKIQQQgjhYCQBFEIIIYRwMJIACiGEEEI4GEkAhRBCCCEcjCSAQgghhBAORhJAIYQQQggHIwmgEEIIIYSDcbJ1AO1Z1SYq+fn5No5ECCGEEI1V9bntyJuhSQLYAtnZ2QCEhITYOBIhhBBCNFV2dja+vr62DsMmJAFsAT8/PwCOHz/uUG+g/Px8QkJCyMrKcqg9FOW85bwdgZy3nLcjyMvLo1u3btWf445IEsAWUKsrh1D6+vo61I1TxcfHR87bgch5OxY5b8fiqOdd9TnuiBz3zIUQQgghHJQkgEIIIYQQDkbz9NNPP23rINozjUbDxIkTcXJyrN50OW85b0cg5y3n7QjkvB3rvKuoFEeeAy2EEEII4YCkC1gIIYQQwsFIAiiEEEII4WAkARRCCCGEcDCOOfKxCRYvXoxWqwXAaDSycOFCqxzT1ixevBiAzMxMABITE+t9vV6vJzExkaioKHQ6HampqYwYMYLo6Girx2opzT2H9n69Z82aRVRUFMOHD68+jyo6na7WY9rj9TYYDOj1epKTk0lNTb2i3F7v9cacN9jfvV7fedvzvV7fedv7vd7Qe9le7/FmU0SdEhISlMTExOrHqampSkxMjMWPaWsWLlxY43FMTIwSGRlZ7zHJycmKTqdTAEWn09X4G7QXzTkHe7jeVed8+U90dHSdx7S3652enq4kJiYqCQkJSkRExBXl9nqvN3Te9nqvN3Te9nqvN3Te9nyvN/Rettd7vCUkAayHVqtVcnNzazzXUM7cnGPaktzcXCUyMrLGOaSnpyuAkpmZWedxycnJV5x3e9Occ2jv11tRrvzHqShKg//k2+v1Tk5OrvWD0d7v9drO2xHu9bqut73f63Wdt73e6415L9v7Pd4cMgawDgaDAaPReEUzOVQ2iVvqmLYoLS0Ng8FQ/biqa8BoNNoqpDbJXq53bGxsjcd6vZ7hw4fbKJrWJ/e63OsNsZfrbc/3en3vZUe+x+sjYwDrcOkb6VJarbbOf47NOaat0Wq15Obm1niu6s1e1xiRKqtWrareWHv79u0kJCRYJ0graso52MP1hprXteqfZWRkZIPH2cP1BrnXLyX3eu3s4XqD/d7rDb2X09LS6jzOnu/xhkgC2ER+fn7k5ORY/Zi2JD4+nsTExFq/CVXR6XRERERU/4PJyclh1qxZJCcnt1aYLWapc2jP1zsuLq7BSQBgH9e7IXKv184err3c6/Z/rzfmveyI9/ilpAu4iZpz4dvzmyUuLo7Y2FhiYmLqfd2l/yAAZs+eTUpKSrv6pmSpc2iv1zsjI6PRr7WH690QuddrZw/XXu51+77XG/tedrR7/HKSANahri4Qo9FYZ1lzjmnLUlJSCAsLa/Amqnrtpaq+ddXVjN4WNfUc7O16JyYmEhYW1qjX2sP1riL3utzrcq/Xrb1d79rey3KP18HWs1DaMq1We8VsuIb+ZM05pi1KTU1VkpOTqx/n5ubWOTMwNzf3ipmDVc+19dljVZp7DvZyvRVFafQSD+35etc3C9ie7/W6zltR7Pter2v2s73f6/Vdb0Wx33u9vveyvd/jzSEtgPVYtGhRjdk+KSkpNb5VGAyG6oUnG3tMe5CRkUFGRgYREREYDAYyMjKIj4+vHgR8+XlrtVoSExNrfCtKSkoiOjq63vEXbUljzsFer3cVg8FQ6zdbe7redXXf2Pu9Xtd52/u9Xtt5O8K93lA3pT3e6w29l+39Hm8OlaIoiq2DaMsuXQU8MzOzxuynpKQkEhISqlcdb8wxbZ3RaCQ0NLTW8R1Vb5XazttoNJKUlFT9ODs7u12dNzR8DvZ4vS8VFhZGcnIyERERNZ63h+ttMBhITExEr9eTkZHBwoULr+gmssd7vb7ztud7vaHrba/3emPe52B/93pj3stgn/d4S0gCKIQQQgjhYKQLWAghhBDCwUgCKIQQQgjhYCQBFEIIIYRwMJIACiGEEEI4GEkAhRBCCCEcjCSAQgghhBAORhJAIYQQQggHIwmgEEIIIYSDkQRQCCGsIDY2ltjY2Fp3J7hcUlISsbGxpKSktEJkQggBTrYOQAghGhIXF4der2f48OGEhYWxfft29Hp9ja299Ho9CQkJREdHVx9XlXy19r6lRqMRrVbb6G2jqs4jLi6uRvxCCGEtkgAKIdo8g8FAenp69ePFixeTkZFRI8FKSUm5ItFbsGABc+bMkaRKCCEuIwmgEKJNy8jIYNGiRTWe2759+xUb2Wu1WnQ6XY3nkpOTrR6fEEK0RzIGUAjRphkMhiuSPb1eT1RU1BWvvTwBFEIIUTtpARRCtGmXd98aDAaMRiORkZE1nr/0sV6vx2AwAH+Mr4PKyRaZmZnExsZW17N9+3bCwsKIiYmpnoSxcuVKli1bdkWX8uLFi9HpdBgMBiIjI69ITBtiNBpJSkpCp9NV1335eQghRGuQBFAI0a5kZGTU2t1bpSrxi4yMJCoqqjoBzMjIIDIyEqPRyKxZs0hOTkan0xEdHU2HDh2AP5LF7du3k5SUxMKFC6vrHTZsWPUxUDlho6kJYFxcHHFxcdVJZFxcnCSAQgibkC5gIUS7kpqaWm/SVNU6p9frayRoOTk56HQ6MjMziYyMrE7kjEYjRqOR2bNn16jj0ta/2NhYhg8fXn2MwWCgY8eOTY49LS2NhIQEDAYDOp1OxigKIWxGEkAhRLui1+sZMWJEneVVyWFiYiKxsbFXPK/X65kzZ07182lpaURGRtZI+DIyMhg+fHj146SkJLRaLUlJSSQlJZGSklKjdbCxEhISWLVqFWFhYQwbNqzJxwshhKVIF7AQot0wGo3VLXz1MRgMtb6u6vhLWwZTU1NrPDYajeTk5FQ/V9Wl3Ng1/eoTGRlJbm4uGRkZzJo1i5SUFFmiRghhE9ICKIRoN9LS0gAaHHuXkpJS3aWr1+trHH95Unj5jOJVq1bVOLZqUeeWCgsLq04mIyIiiI6ObvUFqoUQoookgEKIdqOh8X9Vtm/fTlRUFAaDAT8/vxrHX548Vk0OufQ1VcdC7cmm0Whs0rZtRqOR6OjoGhNXGtOSKYQQ1iJdwEKINm/x4sVkZ2eTlJSEn58fcXFxzJkzp86WwDlz5lQncJd2sRoMhhqLSteWhMXGxpKRkVHj2O+++656CRioXHS6KV23Wq2WESNGVCeNOTk5LFu2rNHHCyGEpakURVFsHYQQQtgTo9FIfHx8k8cNxsXFWWSsoRBCNES6gIUQQgghHIwkgEIIIYQQDkYSQCGEsIKUlBRiY2MxGo0NvjYpKYlZs2a1QlRCCFFJxgAKIYQQQjgYaQEUQgghhHAwkgAKIYQQQjgYSQCFEEIIIRyMJIBCCCGEEA5GEkAhhBBCCAfz/6h5dxtl9r3gAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<IPython.core.display.Image object>"
      ]
     },
     "execution_count": 15,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "from IPython.display import Image\n",
    "Image(\"kundur_full_out_1.png\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Note: `andes plot` does not display the image correctly in iPython or notebook as of now. Solutions are welcome."
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.7.6"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}
